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Abstract 

Elaborated calculations of the shear and the bulk viscosities in the hadron gas, using the ultra- 
relativistic quantum molecular dynamics (UrQMD) model cross sections, are made. These cross 
sections are analyzed and corrected, and the consequences of the corrections for the freeze-out 
temperatures are discussed. A special treatment of the resonances is implemented additionally. 
The previously considered approximation of one constant cross section for all hadrons is justified. 
It's found that the bulk viscosity of the hadron gas is much larger than the bulk viscosity of the 
pion gas while the shear viscosity is found to be less sensitive to the hadronic mass spectrum. The 
maximum of the bulk viscosity of the hadron gas is expected to be approximately in the tem- 
perature range T = 150 — 190 MeV with zero chemical potentials. This range covers the critical 
temperature values found from lattice calculations. A misleading viewpoint, resulting in mistakes 
in the calculations of the bulk viscosity in some papers, is commented on. Doing this, a generalized 
Chapman-Enskog procedure, taking into account deviations from the chemical equilibrium, is out- 
lined. Some general properties, features, the physical meaning of the bulk viscosity and some other 
comments on the deviations from the chemical equilibrium accomplish this discussion. Analytical 
closed-form expressions for the transport coefficients and some related quantities within a quite 
large class of cross sections can be obtained. Some examples are explicitly considered. Comparisons 
with some previous calculations of the viscosities in the hadron gas and the pion gas are done. 
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I. INTRODUCTION 



The bulk and the shear viscosity coefficients are transport coefficients which enter in the 
hydrodynamic equations, and thus are important for studying of nonequilibrium evolution 
of any thermodynamic system. 

There are two more additional reasons to study the shear viscosity. The first one is 
the experimentally observed minimum of the ratio of the shear viscosity to the entropy 
density rj/s near the liquid-gas phase transition for different substances, which may help 
in studying of the quantum chromodynamics phase diagram and finding of the location of 



the critical point 
models, see e. g. 



bound 2 77 js 
ones see [121 




2J 1 . Such a minimum was observed in theoretical results in several 
5(. For a counterexample see [6[ and references therein. The second 
reason is the calculation of the i]/s in strongly interacting systems, preferably real ones, 
to compare physical inputs which provide small values of the rj/s. The conjectured lowest 
= [TP ] was violated with different counterexamples. For some reasonable 
13. Also see the recent review [14|. The bulk viscosity, being very sensitive 
to violation of the equation of state and being connected with fluctuations through the 
fluctuation-dissipation theorem [15J], can have a maximum near a phase transition 

Bis- 

18]. In |i| and [18| sharp maxima were observed in the bulk viscosity £ and the ratio £/s 
in the linear cr-model for the vacuum a mass 900 MeV. Decreasing the vacuum a mass 
the maximum eventually disappears. Any maximum of the £/s was not observed in the 
large-N limit of the linear a- model in the [18|] 3 . Any maximum of the £ was not observed 
in the large-N limit of the 1 + 1-dimensional Gross-Neveu model [2(| (see also Sec. [IT] for 



1 Fireballs, created in heavy ion collisions, have finite sizes and finite times of existence of their thermal- 
ized part. This puts important restrictions on detection of the critical fluctuations of thermodynamic 
functions [3]. Because of this it's also important to consider nonequilibrium dissipative corrections and 
nonequilibrium phenomenons like critical slow down/speed up. 

2 In jj the bound coming from the Heisenberg uncertainty principle was obtained for the rj/s. However, it 
was obtained using a formula, which is justified in rarified gases with short-range interactions. It's well 
known already from the nonrelativistic kinetic theory that dense gases get corrections over the particle 
number densities (see e. g. Section 18), corresponding to more than binary collisions, and in very dense 
gasses this bound can be quite inaccurate. In liquids and other substances the mechanism of appearance 
of the shear viscosity may be different (see [9] for a review). In particular, the shear viscosity of water 
can be very well described by a phenomenological formula with an exponential dependence on the inverse 
temperature, see e. g. fl~C)j |. 

3 There in the case of crossover (in the calculations with massive pions at T = 0) the squared speed of 
sound without thermal corrections to the mass coincides with the analogical squared speed of sound in 
the charge-neutral (or chargeless) system, namely, the pion gas squared speed of sound (see fig. [3] for 
comparison). Thermal corrections to the mass make the squared speed of sound to cross the value 1/3, 
and the bulk viscosity become equal to zero approximately at the same temperature. Note that the pion 
gas bulk viscosity has a maximum at T = 40 — 50 MeV, defined by the pions' masses, but it's not related 
to any phase transition. In the hadron gas with a reach mass spectrum the squared speed of sound has 
no such a behavior before and after the introduction of the interactions, see UJ1. Also, as is explained 
in Sec. UH at low enough temperatures the bulk viscosity source term should be changed on that one of 
a system with only the elastic collisions (discarding the inelastic collision integrals for self-consistency of 
the Boltzmann equation). 
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comments). 

Whether one uses the Kubo 4 formula or the Boltzmann equation one faces nearly the 
same integral equation for the transport coefficients j23l-[25|. The preferable way to solve it 
is the variational (or Ritz) method. Due to its complexity the relaxation time approximation 
is used often in the framework of the Boltzmann equation. Though this approximation is 
inaccurate, does not allow to control precision of approximation and can potentially lead 
to large deviations. The main difficulty in the variational method is in calculation of colli- 
sion integrals. To calculate any transport coefficient in the lowest order approximation in 
a mixture with a very large number of components N' (like in the hadron gas) one would 
need to calculate roughly iV' 2 12-dimensional integrals if only the elastic collisions are con- 
sidered. Fortunately, it's possible to simplify these integrals considerably and perform these 
calculations in a reasonable time. 

This paper contains calculations of the shear and the bulk viscosity coefficients for the 
hadron gas using the (corrected, see Sec. IIIII) UrQMD cross sections. The calculations are 
done in the framework of the Boltzmann equation with the classical Maxwell-Boltzmann 
statistics, without medium effects and with the ideal gas equation of state. The Maxwell- 
Boltzmann statistics approximation allows one to obtain some relatively simple analytical 
closed-form expressions. Originally the calculations in the same approximations for the 
hadron gas but with one constant cross section for all hadrons were done in (26[. The 
deviations in the worst cases are relatively small. In that paper some analytical formulas 
of the viscosities for 1-, 2- (explicitly) and N-component (up to solution of the matrix 
equation) gases with constant cross sections were obtained. Analogical formulas can be 
written down for quite a large class of non-constant cross sections, in particular, for the ones 
which appear in the chiral perturbation theory. The final expressions may become somewhat 
more cumbersome; anyway this is better than numerical integration at least in the speed 
of the computation. Explicit formulas for the viscosities with the elastic pion-pion isospin 
averaged cross section and somewhat more general one are obtained in the present paper. 
The results of the j26| are partially reproduced in the present paper, improving the text and 
adding more detailed explanations. The presented calculations can be considered as quite 
precise ones at low temperatures where the elastic collisions dominate and the equation of 
state is close to the ideal gas equation of state. At higher temperatures the calculations with 
the total cross section are expected to give the qualitative description. 

For comparison the calculations of the viscosities are performed for the pion gas (through- 
out the paper the chemical potentials are equal to zero if else is not stated). The results 



are relatively close to the results in [27j, [28|. There the calculations are made in the same 
approximations except for the [28] , where the Bose- Einstein statistics is used instead of the 
Maxwell-Boltzmann one. The discrepancies from the used classical statistics are not large 
at zero chemical potential and become larger as the chemical potential grows (see Sec. IIVI 
for the errors and comments). The comparison is made with the results of the 27], see 



fig. [TJ The discrepancies up to a factor of 3.5 for the bulk viscosity and up to a factor of 
2.5 for the shear viscosity come most probably from somewhat different tttt elastic plus the 
quasielastic (through the intermediate p-resonance) cross section of the |27j 5 (the averaging 



21 



22|. 



4 The Kubo formulas are distinguished from the Greet-Kubo formulas e. g. in 

5 The author could not reproduce this plotted total cross section by its formula. In fact it was approximately 
2.6 times larger. But the plotted total cross section is quite close to the isospin averaged (corrected) 
UrQMD 7T7T total cross section. A notable deviation is only at y/s < 0.5 GeV, when the UrQMD cross 
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over the scattering angle is expected to give small errors; see also comments below the for- 
mula (17.25P ). The minima of the shear viscosities near T = 60 — 70 MeV are attributed to 
the peaks from the p-resonances in the tttt cross sections. It's not noticeable in the figure 
for the dash-dotted line. Nonzero values of the bulk viscosities and theirs maxima are solely 
due to the masses of the pions. The paper 27j implements also the isospin averaged current 
algebra elastic cross sections. These cross sections can be reproduced in the lowest order 
in the chiral perturbation theory (29[. They obviously have quite large deviations from the 
experimental data at high enough energies and wrong y/s asymptotic dependence, which can 
be seen from the comparison of them with the isospin averaged elastic plus the quasielastic 
experimental cross sections in the 27 
constant 5 mb 130, 31 



The elastic tttt cross sections are rather close to the 
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FIG. 1: The viscosities of the pion gas as functions of the temperature at zero chemical potentials. 
The dotted and the dash-dotted lines designate correspondingly the bulk viscosity time 10 3 and the 
shear viscosity, calculated in the 3-rd order in the 271 ] ■ The solid and the dashed lines designate 
correspondingly the bulk viscosity time 10 3 and the shear viscosity, calculated at the same orders 
using the (corrected, see Sec. IIIip UrQMD cross sections. 



In several papers the bulk viscosity £ was calculated for the pion gas, using the chiral per- 
turbation theory (or the unitarized chiral perturbation theory) and some other approaches, 



with quite large discrepancies between the quantitative results. In [32| the calculations were 
done by the Kubo formula in a rough approximation. There the number- changing 2 -H- 4 
processes were neglected too, and the non-vanishing value of the bulk viscosity is obtained 
due to a trace anomaly and the pions' masses. At small temperatures, where the effects 
of the trace anomaly are small, the magnitude of the bulk viscosity is large in compare 
to the results of this paper and the [27|, |28|. For example, at T = 20 MeV (60 MeV) it's 
larger approximately in 39 (8) times than the bulk viscosity in this paper. The maximal 
values differ in approximately 24 times. In 33] the calculations are done in the framework 



of the Boltzmann equation and have a divergent dependence of the £ for T — > because 



section becomes 1.5 — 2 times smaller. At 1.2 GeV < y/s < 1.9 GeV the UrQMD cross section is a little 
larger instead. 
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of remained weak 2 f-> 4 number- changing processes (at T = 140 MeV the bulk viscosity is 
nearly 57 times larger than the bulk viscosity calculated in this paper). This dependence 
should change at low enough temperatures, T = 180 MeV or higher ones for the pion gas, 
see Sec. |TT] for explanations. Joining the results of the calculations at low and high tem- 
peratures, the function £ (T) may turn out to be not continuous at the middle temperatures 
(which is not a physical effect, see Sec. [TTJ) , and the smooth function £(T) is to be obtained 
through some interpolation. In [34j the bulk viscosity was calculated in the framework of the 
Boltzmann equation with the ideal gas equation of state and only the elastic collisions taken 
into account. The Inverse Amplitude Method was used to get the scattering amplitudes of 
the pions. The quantitative results are close to the results in this paper (discrepancies up 
to a factor of 2.7). In [35] the calculations are done in the framework of the Boltzmann 
equation for the massless pions. There the bulk viscosity increases rapidly so that the ratio 
£/s increases with the temperature. 

Calculations of the shear viscosity in the hadron gas with a large number of components 



were done in [36j , using some approximate phenomeno logical formula, and in [19| , using the 
relaxation time approximation. These results are in good agreement with the calculations 
of this paper. Hence, as long as the ratio £/r] calculated in the [l9| for the free massive pion 
gas is 8 — 58 times larger (in the temperature range 0.02 GeV < T < 0.14 GeV with the 
deviations growing as the temperature decreases) than the one calculated in this paper, one 
can suspect that the difference comes from the bulk viscosity because of the used relaxation 
time approximation 6 and likely not conserved particle numbers at low temperatures, pro- 
vided that the SHMC model's cross sections, used in the [ll, don't have large deviations 
from the UrQMD cross sections or the experimental data, which seems to be the case. Also 
note that the results in the [l9| for the free particles and the SHMC model don't differ very 
much. These facts may explain why the £/s of the hadron gas in the [l9| is 3.65 — 11.3 times 
larger (in the temperature range 0.1 GeV < T < 0.18 GeV) than the £/s calculated in this 
paper. At the low temperature T = 100 MeV and the vanishing chemical potentials it is 
11.3 times more (at the same temperature the factor is 8.2 for the case of the pion gas). In 



37J the calculation of the bulk viscosity is done for the hadron gas (with an excluded- volume 
equation of state) with the masses less than 2 GeV using some special formula, obtained 
though some ansatz T^|. Its quantitative accuracy has not been clarified. The ratio £/s in 



the [37J deviates from the £/s of this paper up to a factor of 1.8 in the temperature range 
0.14 GeV < T < 0.18 GeV and is different on 4% at T = 0.14 GeV. 

Also the shear viscosity has been calculated using the Kubo formula (or the Green-Kubo 



formula) in a gas of mesons and their resonances [21] . There the UrQMD simulations are 



performed to calculate the energy-momentum tensor, used in the calculations by the Kubo 



formula. The rj in the [21[ is 1.14 — 1.77 times smaller then the rj for the hadron gas in this 
paper. At T = 0.15 GeV it is 1.77 times smaller. In [38[ similar calculations, using the 
URASiMA event generator, are done for the shear viscosity with close results. 



6 In the relaxation time approximation the bulk viscosity source term is treated somewhat differently: the 
£ becomes proportional to the integral of the squared source term (times some functions of momentum) 
and not to the square of the integrated source term (times some functions of momentum) . Note that in 
the [32I the used formula has this relaxation time approximation form. Also there the source term is the 
one of a system with the inelastic processes. These facts could help to understand the enlarged values of 
the bulk viscosity. Not small quantitative discrepancies can be noticed between the calculations of the 



and the 181 
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The structure of the paper is the following. A misleading viewpoint on the bulk vis- 
cosity, connected with the inelastic processes, is commented on in Sec. [Til together with 
some properties, features and physical meaning of the bulk viscosity. In that section some 
questions concerning the deviations from the chemical equilibrium are addressed too. Sec. 
II III contains some comments on the constant cross sections, which are used in approximat- 
ing calculations, and some other general comments on cross sections. Also it contains a 
description of the UrQMD cross sections, which are used in the main calculations, together 
with their analysis, corrections and the consequences of the corrections for the freeze-out 
temperatures. The applicability of the used through the paper approximations is discussed 
in Sec. IIVI The system of the Boltzmann equations, its solution and formal expressions of 
the transport coefficients can be found in Sec. |Vj The numerical calculations for the hadron 
gas are presented in Sec. |VTJ In Sec. IVII Al analytical results for the single-component gas 
are presented. In particular, an analytical expression for the first order single-component 
shear viscosity coefficient with constant cross section, found before in j39|, is corrected while 
the bulk viscosity coefficient remains the same. The nonequilibrium distribution function 
in the same approximation is written down. Also the viscosities with some non-constant 
cross sections are written down. Some analytical results for the binary mixture with con- 
stant cross sections are considered in Sec. IVII Bl Integrals of source terms needed for the 
calculation of the transport coefficients can be found in Appendix [A] The general entropy 
density formula can be found in Appendix [Bj It is used in the numerical calculations for the 
hadron gas. Transformations of collision brackets, being the 12-dimensional integrals which 
enter in the viscosities, and some analytical formulas for them can be found in Appendix 
O The closed-form expressions for collision rates, mean free paths and mean free times are 
included in Appendix [D] 



II. SOME FEATURES AND PROPERTIES OF THE BULK VISCOSITY 

First, it should be reminded that the transport coefficients are defined as coefficients next 
to their gradients in the formal expansion of the energy-momentum tensor and the charge 
density flows over the gradients of the thermodynamic functions or the flow velocity (see e. 
g. [io[, Section 136). The Kubo formulas are not definitions of the transport coefficients, as 
one might think. They may introduce some assumptions. In particular, the Kubo formulas 
in the form as in the [23| have zero frequency and zero momentum limits, which neglect 
finite size and finite time effects. Zero momentum limit implies the thermo dynamical limit. 
This limit is needed to avoid possible nonphysical contributions from inappropriate choice 



of a current and an ensemble [22]. The Kubo formulas in the form as in [38|, |4l| suppose 
thermal equilibrium in the initial moment of time t = — oo. So that any infinite space-time 
scale cannot be connected with the transport coefficients by their definitions. 

The Boltzmann equations will be used in what follows. In the case of the elastic collisions 
they can be derived from the Liouville equation 7 in the approximations n&r^ <C 1 (r^ 
is the effective radius of two particle interactions between the particles of the species k 
and /) that is for rarified gases with short-range interactions 8 . Also the linear integral 



7 The Boltzmann equations can also be derived for the case of the inelastic collisions from some physical 



considerations, see [42j (Chap. I, Sec. 2). 



8 This is the case of interest. Coulomb interactions can be neglected in heavy ion collisions at all the 
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equations for the viscosities and other transport coefficients, derivable from the Boltzmann 
equation, can be obtained (with some corrections) from the perturbative calculations for 
quantum field theories at finite temperature (including the inelastic processes) using the 
Kubo formulas [23|, |43 , which justifies application of the Boltzmann equation when the 



inelastic processes are present. 

The bulk viscosity can reveal itself only when there is a nonzero divergence of the flow 
velocity. This nonequilibrium perturbation should not be confused with another possible 
independent perturbation (as was done in several papers, some of which are mentioned below; 
the roots of the misleading viewpoint, perhaps, can be found in [40|], Section 81). Namely, 
this is the homogeneous perturbation. It can be both the chemical and the kinetic one 9 . 
Then it can be generalized and made dependent on the coordinate. It just should not be 
proportional to any gradient. Then the constraints of the local conservation laws should be 
imposed on these perturbations. The perturbations for the particle numbers should be such 
that don't violate conservation of all charges. Considering the case of homogeneous chemical 
perturbation in a gas with fixed volume, one concludes that the temperature should change 
with time, being some energy per particle. So that energy conservation should be obtained 
varying the temperature. Adding an infinitesimal correction to the temperature one gets 
a perturbation of the form Cp^U^. Such perturbations don't contribute to all collision 
integrals. To describe purely chemical perturbations they have to be chosen in the form of 
the momentum- independent terms (except for the Cp^U^ terms), otherwise there will be 
contributions from the elastic collision integrals. Such perturbations can be considered as 
chemical potentials-like ones (being small, one can expand the distribution functions over 
them and get these momentum-independent terms) with the arguments for maximization of 
the entropy. To find the evolution of these terms they should be separated. Let's write this 
in some formulas. Multicomponent gas distribution functions with the leading perturbations 
can be represented in the form (detailed definitions can be found in Sec. IV Aj) 

fk = /i 0) (l + + <Pk) ~ /f (1 + 0k + <Pk), \0k\ < 1, \<Pk\ < 1, (2-1) 

where (fk are the perturbation due to the gradients and <pk are the chemical perturbations 10 . 
Following steps of Sec. IV A[ one can get the following linearized equations from the Boltz- 
mann equations: 

Pt(U,D + V M )/f + fi°V k (U,D + V M )ft « -/PAN - /f 4" e 'N, (2-2) 

where and D^ el is the sum of the linearized collision integrals (divided on the —f^) of 
all the processes and of the inelastic processes correspondingly. The 2-nd order gradients 
and the squared 1-st order gradients are neglected in the 1. h. s. because they are of the 



considered energies in this paper. 
9 The inclusion of this kinetic perturbation is similar to the inclusion of the chemical one so that it is omitted 
for simplicity below. Usually this perturbation should fade first because also the inelastic processes are 
responsible for the relaxation of the momentum spectra. However, see comments for 1 + 1-dimcnsional 
systems below. 

Note that if the k-th species have conserved particle numbers, then the nonequilibrium chemical potential 
is nonphysical or redefining the usual (thermodynamic) chemical potential. 
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next order 11 and should be cancelled in the next iteration by the next corrections to the 
distribution functions. Also the smallness of the (fk is used. If the spatial covariant gradients 
V At <^fe(t = 0) (at the initial moment of time) are of the same order as the gradients of the 
thermodynamic functions or the flow velocity, then the V^ipk terms in the 1. h. s. should 
be retained 12 . The covariant temporal derivatives Dipk are needed to describe the temporal 
evolution of the (p>k- Then the equations ( 12. 2 p can be split onto the separate equations for 
the ifik and the ipk 

p£U,D^ k + P y^ h » -£rW (2.3) 



Pp.DfW + rfyj* „ -f^C k [<p k }. (2.4) 

The equations ( 12. 2 p can be split within the framework of the perturbation theory over the 
gradients. Let's consider also the condition (p k <C (p k in the (12. ip . Then neglecting the ip k i n 
the ( 12. ip and repeating the steps of Sec. IV A\ one can get the following linearized equations: 

flU^fa a (2.5) 

The equations ( 12. 5 p are precise in the homogeneous case (the approximation is only from 
the linearization). The 1-st order gradients and the V M £>& are neglected. Then using the 
and the ([2~2]) . one can get 

(P£U»D + ^V,)/i 0) + fPrtVtfu « -A 0) Ck[<Pk). (2-6) 

Solving the system of equations (12. 5p in the local rest frame, one gets the leading exponen- 
tial fading dependencies on time 13 (in a covariant form this should be an explicit space-time 



dependence). Such dependencies were obtained in some previous studies, see e. g. [45l. 146 
The equations (I2.6P are different from the ones obtained from the common Chapman-Enskog 
procedure (see e. g. (42[, Chap. V) because of the V^fik terms. The contributions from 
the small chemical perturbations can be neglected in the considered order in the transport 
coefficients because they are multiplied on the 1-st order gradients. The V '^(pk terms can 
be cancelled, introducing terms proportional to the V ^(pkit = 0) into the ipk terms. If the 
spatial distributions of the (p>k{t = 0) are such that V^kit = 0) are of the 2-nd or a higher 
order, then the V^k can be neglected. This assumption or approximation is used in the 



11 The question of validity of this expansion over the gradients (which coincides with the usual order counting 
in the formal expansion over the gradients in the hydrodynamics) for some profiles is not discussed in this 
paper. 

It's a reasonable assumption in the case when the hydrodynamical description is applicable. For example, 
the chemical perturbations can be a result of a fast previous expansion (faster than the chemical equilibra- 
tion). Then the inhomogeneities of the chemical perturbations should correlate with the inhomogeneities 
of the thermodynamic functions, the flow velocity or it's divergence. 
13 If the expansion rate is much larger than the collision rates of the inelastic processes (e. g. because of 
a substantial decrease of the temperature), then the chemical perturbations should enlarge instead. If 
the r. h. s. of the (|2 .3[) is smaller than the second term of the 1. h. s., then one can consider another 
approximation, when the fc-th species particle numbers are conserved. Then the chemical perturbation 
becomes an addition to the thermodynamic chemical potential. 
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calculations of this paper. In the linear response theory one can also introduce indepen- 
dent small chemical perturbations with the same conclusions for the 1-st order transport 
coefficients and find evolution of the perturbations with time. 

Note that the deviation from the chemical equilibrium itself is not necessarily a source 
of the bulk viscosity, as is stated in [47]. If the bulk viscosity is not equal to zero only 
because of the particles' masses and they are tended to zero, the bulk viscosity source term 
and the bulk visc osity tend to zero even if there are inelastic processes (see the end of 
Sec. IV All . In the [47| the independent chemical perturbations and the perturbations due 



to the gradients were just connected through the perturbations of particle numbers, and 
the bulk viscosity became proportional to the chemical relaxation time. Formally infinite 
chemical relaxation time doesn't imply any divergencies in the chemical perturbations 
but rather approximation of conserved particle numbers. Note that the dependence on 
the strength of the inelastic processes is different for the chemical perturbations and the 
perturbations due to the gradients. Increasing the strength of the inelastic processes the 
chemical relaxation time decreases. And the gradients' relaxation time increases, because 
the transport coefficients, at least in rarified gases with short-range interactions, roughly 
speaking, are inversely proportional to the integrated cross sections (in an ideal liquid the 
gradients' relaxation time is infinite). What happens with the bulk viscosity if the inelastic 
processes become weaker is discussed below. 

Making the inelastic processes weaker in compare to the elastic ones, the bulk viscosity 
eventually gets a formal dominant contribution from them because of the approximate zero 
mode(s) [23j, connected with possible conservation of particle number(s) 14 . As long as 
it's clear that the bulk viscosity is not responsible for the chemical equilibration, it's also 
clear that there may be the approximation of conserved particle numbers if the momentum 
spectrum, as well as the gradients, can relax by means of only the elastic collisions (which 
is usually the case) and the elastic processes make a dominant contribution to the collision 
rates. The question is only at what concrete temperature does this approximation sets in. 
Let's make an illustrative example (resembling the case of gases at room temperatures) of 
what nonphysical contributions one can get from formally remained weak inelastic processes. 
Consider infinitely weak inelastic processes and the perturbation of the flow velocity such 
that the energy-momentum tensor gets a sizable contribution from the bulk viscosity term, 
not large in compare to the pressure (cf. (I5.9p . (15.161) ) to remain the perturbation theory 
applicable. Then it's obvious that this contribution is not physical because it is created 
by the practically absent processes and the infinitesimal perturbation of the flow velocity. 
Instead, this system is practically described by the equilibrium thermodynamic functions. 
This also answers positively the question whether the thermodynamic chemical potential can 



14 



If the particles involved into the inelastic processes are massive, then the formal dominant contribution 
is the exponential one over the temperature and grows as the temperature decreases. If the particles are 



massless or approximately massless, as in high-temperature QCD 48[, then a more complicated situation 



can occur, and one may need to compare some differences of processes' rates (and not just equilibrium 



collision rates), arising in the collision matrix (C^ n in assignments of the [3|)- Under the same pair 
of used test-functions (indexed by to, n in the Cf m ), and for the same pair of particle species, smaller 
differences of processes' rates can be neglected. Comparing among different pairs of test-functions the 
smallest nonzero contributions dominate, or rather as can be obtained directly from the inverted collision 
matrix. 
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be introduced for approximately conserved particle number in principle. As far as the author 
knows, the first correct comment (albeit somewhat inaccurate) on this issue can be found in 



the [24j. However, note that in fact there is no divergent mean free paths, corresponding to 
the inelastic processes in this case. They are cut by the mean free paths, corresponding to 
the elastic processes (and the overall collision rate have the dominant contribution from the 
elastic collisions). So that it may be not necessary for the chemical relaxation time to be 
much larger than any relevant time scale (like the gradients relaxation time or the time of 
existence of the thermal part of the system) to switch off the inelastic processes. That's why 
a criterion based on comparison of collision rates of elastic and inelastic processes can be 
considered to switch off the inelastic processes. Such a comparison is done in the UrQMD 
studies of the hadron gas in [49[ (see Sec. [IV] for farther discussions). According to [50], 
the chemical relaxation time of the 2 -h- 4 processes in the pion gas is much larger than the 
thermal relaxation time. And e. g. at T = 180 MeV the chemical relaxation time is equal 
to 40 fm, which is larger than the typical lifetime of the thermal part of the expanding 
fireball (see e. g. [49[). So that it's the inelastic 2 <H- 4 processes which should be neglected 
in the pion gas at T = 180 MeV or even higher temperatures, which wasn't done in the 
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The bulk viscosity source terms increases substantially if particle numbers are not con- 
served (cf. (I7.17p . (j7.18P ; in mixtures these particle numbers should also be not small). This 
reflects additional fluctuations from not conserved particle numbers. Though the inelastic 
processes have to be effective enough to consider the approximation of not conserved particle 
numbers. Perhaps, the point at which the bulk viscosities in the different approximations 
cross can provide a criterion for switching on/off the inelastic processes. If this is not so, 
then one would have to make some interpolation in the intermediate region 15 . Note that e. 
g. in the calculations by the Kubo formulas through the direct calculations of the energy- 
momentum tensor as in the 21] it's not needed to use the approximation of conserved or 
not conserved particle numbers (which defines the number of independent thermodynamic 
chemical potentials, through which the chemical potentials of all particles are expressed, cf. 
( 15. 7p ). There the energy- momentum tensor should be a smooth function of time and the 
thermodynamic functions as long as the inelastic processes fade smoothly. Then the bulk 
viscosity should be a smooth function of the temperature and particles' chemical potentials 
regardless of the number of the independent chemical potentials. 



In the [48[ a bottleneck for the relaxation to equilibrium characterized by the bulk vis- 
cosity due to the weakest processes' rates is assumed. Instead, there are rather dominant 
contributions from some test-functions 16 (as is commented in the footnote [14")) . which should 
not be specially treated though, except for the ones which are the approximate zero modes 
making a dominant contribution. A similar dominance 17 is present also in other transport 
coefficients, in particular, when there is only one type of processes. Although in QCD at 
high enough temperatures the equilibrium 2 <H- 2 elastic collisions rate is parametrically the 



15 Perhaps, the bulk viscosity calculated without constant test-functions (except for zero modes of the 
inelastic collision integrals, used to conserve charges) can provide a good interpolation. 

16 Not a bottleneck from some perturbations, because one actually doesn't have a choice in the form of 
the momentum dependence of the perturbations corresponding to the transport coefficients. The kinetic 
perturbation can be of different forms of the momentum dependence. 

17 Another similar dominance can exist from particle species interacting weakly with all particles. 
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largest one 18 , 0(a s T), because of cancellations the momentum transfer takes place with the 
rate 0{pt 2 s \n{l/a, s )T), which is parametrically smaller than the particle number change rate 
0{o?J 2 T) from the effective "1 -H- 2" processes. This provides an example when the equilib- 
rium collision rates may differ substantially from the relevant collision rates. The " 1 -<->■ 2" 
processes provide small chemical relaxation time in compare to the thermal relaxation time, 
which justifies the approximation of not conserved particle numbers and the enhancement 
of the bulk viscosity from the source terms at least at small enough a s , whereas the con- 
tributions to the bulk viscosity from the collision integrals of the "1 -H- 2" processes are 
suppressed at small enough a s (the inelastic 2 -h- 2 processes are not suppressed, but they 
are of the order 0(a 2 hi(l/a s )T)). To avoid misunderstanding it may be mentioned that 
taking the total collision rate of the "1 ■<-»■ 2" processes as formally infinite by taking the 
corresponding matrix elements as formally infinite ones, one gets zero bulk viscosity and 
zero mean free paths as long as both the gluons and quarks take part in these processes (see 
also footnote . 

In the case of a 1 + 1-dimensional single-component gas the elastic collisions cannot 
result in the relaxation of the momentum spectra and, hence, cannot stimulate the system 
to evolute towards equilibrium 19 . As a result, the exponentially divergent bulk viscosity 
was obtained in the paper |20|. Considering again the example about the infinitely small 
perturbation of the flow velocity and assuming also a finite size of the system, it's again 
obvious that the weak inelastic processes may make nonphysical contributions (in this case 
the mean free path is formally cut by the system's size). If this is the case, then the 
hydrodynamical description becomes inapplicable, and might use simulations of particles' 
collisions or the Boltzmann equations in the approximation without collisions (on a time 
scale much smaller than the chemical relaxation time). If the 1 + 1-dimensional description 
is only an approximate one (that is with small angle elastic scatterings in higher dimensions), 
the relaxation of the momentum spectrum by the elastic collisions should be considered. And 
if a 1 + 1-dimensional gas has at least two components with different masses, then a nontrivial 
momentum exchange in the elastic collisions is possible. This results in the possibility of 



the relaxation of the momentum spectra by only the elastic collisions 52 . 

Let's summarize this section with formulation of the physical meaning of the bulk viscos- 
ity. The bulk viscosity reflects deviation of the value of the pressure from its local equilibrium 
value (as can be seen from the ( I5.16P ). appearing when the system expands/compresses, 
because of the delay in the equilibration. The bulk viscosity is not responsible for the 
restoration of the chemical or the kinetic equilibria - it's responsible for the relaxation of the 
divergence of the flow velocity. If there are inelastic processes, then the particle numbers also 
get nonequilibrium contributions (cf. (|5.3|) . (I5.26p . (I5.57P ) such that the charge is conserved 
locally (cf. (I5.63P ) 20 . Though these contributions together with the contribution to the pres- 
sure may become nonphysical because of the approximate zero modes (if such ones appear 
in the calculations). The magnitude of the bulk viscosity changes from theory to theory. 
Under some quite general assumptions a nonzero value of the bulk viscosity can be connected 



18 The estimate can be easily inferred from 5l| . 

19 There are forward scatterings and momentum interchange. As long as the particles are not distinguishable 
the momentum interchange from the elastic collisions is equivalent to the forward scatterings or absence 
of the elastic collisions at all. 

20 One should keep in mind that while studying the chemical perturbations (pu through the thermodynamic 
functions first the contributions from the transport coefficients' terms should be subtracted. 
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with violation of the scale invariance due to a nonzero value of the energy-momentum tensor 
53l . l54j . Of course, the beta function can contribute to the energy-momentum tensor and 
the bulk viscosity too 23 . 



III. THE HARD CORE INTERACTION MODEL AND THE URQMD CROSS 
SECTIONS 



In a non-relativistic classical theory of particle interactions there is a widespread model, 
used in approximate calculations, called the hard core repulsion model or the model of hard 
spheres with some radius r. For its applications to the high-energy nuclear collisions see 



36[ and references therein. The differential scattering cross section for this model can be 
inferred from the problem of scattering of point particle on the spherical potential U(r) = oo 
if r < a and U(r) = if r > a [jjBj]. In this model the differential cross section is equal to 
a 2 /4. To apply this result to the gas of hard spheres with the radius r one can notice that 
the scattering of any two spheres can be considered as the scattering of the point particle 
on the sphere of the radius 2r, so that one should take a = 2r. The total cross section otot 
is obtained after integration over the angles of the r 2 dfl which results in the Otot = 47rr 2 . 
For collisions of hard spheres of different radiuses one should take a = + ri or replace the 
r on the !*±&: 

<?tot,ki = ir{r k + ri) 2 . (3.1) 

The relativistic generalization of this model is the constant (not dependent on the scattering 
energy and angle) differential cross sections model. 

The hard spheres model is classical, and connection of its cross sections to cross sections, 
calculated in any quantum theory, is needed. For particles, having a spin, the differential 
cross sections averaged over the initial spin states and summed over the final ones will be 
used. 21 If colliding particles are identical and their differential cross section is integrated 
over the momentums (or the spatial angle to get the total cross section) then it should be 
multiplied on the factor | to cancel double counting of the momentum states. These factors 
are exactly the factors next to the collision integrals in the Boltzmann equations (15.191) . 
The differential cross sections times these factors will be called the classical differential cross 
sections. 

The UrQMD cross sections are used in the numerical calculations of Sec. IVIP 2 . These 



cross sections are described in |30l. l3ll|. More details can be found in the UrQMD program 



It's assumed that particle numbers of the same species but with different spin states are equal. If this were 
not so then in approximation, in which the spin interactions are neglected and probabilities to have certain 
spin states are equal, the numbers of the particles with different spin states would be approximately equal 
in the mean free time. With equal particle numbers their distribution functions are equal too. This allows 

one to use the summed over the final states cross sections in the Boltzmann equations. 
Very high energy dependence of any used UrQMD cross section is not important because of the exponential 
suppression e _v ^/ T . The used cross sections were cut on the ^/s = 5 GeV and were continued by a 
corresponding constant continuously at higher energies. At small enough momentums there is another 
somewhat weaker suppression. The momentum space density of each particle provides p 2 suppression. 
This may (partially) suppress some deviations from the experimental data of some UrQMD cross sections 
(like for the Ap pair) at \fs ~ rrik + mi- To estimate at what temperatures some discrepancies in cross 
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codes. Below there is some description mainly of what is different or new. 

The UrQMD cross sections are averaged over the initial spin states and summed over 
the final ones. As long as the UrQMD cross sections are total ones (integrated over the 
scattering angle), the factors 7^ are already absorbed into them (in what follows only such 
cross sections will be considered in this section tacitly). Dividing them on the 4ir, one gets 
the classical differential cross sections, averaged over the scattering angle. 

The UrQMD codes (version 1.3) were modified to get accurately tabulated (with a step of 
25 MeV) cross sections. Resonances' masses and widths (they are tuned in their uncertainty 
regions to describe the experimental data better), used in the UrQMD codes, have somewhat 
different values than the ones in the [30]. Influence of variation of these parameters was 
studied in (56[. The UrQMD codes implement somewhat different averaging of the c. m. 



momentums over the resonances' masses 23 than in the papers 1301, 31]. It was found that 
using the resonance dominating cross sections from the papers [30, |3l| some of these cross 
sections could have a large rise at small c. m. momentums if constant widths are used in 
the calculations of the averaged c. m. momentums in the energy dependent widths. So 
that one should be aware of this fact 24 . The UrQMD codes have a low energy cut-off at 
^/i ~ rrik + mi + 0.01 GeV (and a similar one over the c. m. momentum if triggered) for 
the resonance dominating cross sections, and no large low energy rise was found there. 

An important ingredient of the UrQMD model is the Additive Quark Model (AQM), 
which is used for unknown cross sections. Universality of hadrons, based on jet quenching 
arguments, is used to support this model. This model describes the experimentally known 
cross sections well at sufficiently high energies. Application of this model is better than 
elimination of the corresponding hadrons, which is the same as equating their all cross sec- 
tions to zero and, hence, exclusion of their contributions from the thermodynamic functions 
(infinite mean paths, no thermalization). 

At this point an interruption should be made to consider some important questions 
related to different types of the UrQMD cross sections. These different types are used due 
to several reasons and are the following: the elastic cross section(s) (ECS(s)), the elastic 
plus the quasielastic cross section(s) (EQCS(s)), the total cross section(s) (TCS(s)) and 
the previous two types with enhanced in some way resonances' cross sections (index "2" is 
appended in the abbreviations). 

Of course, the system of the Boltzmann equations would have a solution with any of these 
cross sections. Usage of the ECSs is completely self-consistent as long as only the elastic 
2 <H- 2 collision integrals are used in the calculations of the viscosities. However, there are 
reasons to consider also the EQCSs. Exactly this type of cross sections, being averaged over 
the isospin, is implemented in [27]. The quasielastic cross sections 2 — > 1 — > 2 can be used 
as rightful contributions to the ECSs in the approximation that the 4-momentum of the 
intermediate resonance does not change (the effects of the exclusion of the resonances as 
independent particles are considered in Sec. I VI p . The mean free paths of the intermediate 
resonances without contributions of the decays, being not equal to zero, also introduce some 
errors, which are neglected. The EQCSs conserve particle numbers, which is consistent 



sections can appear one can equate the -^i to the sum of the averaged one-particle energies eu (|5 . 1 2[) of 
the two colliding particles. 

Averaged powers of the momentums are used, not powers of the averaged momentums. 

It may be mentioned that one should be also aware of possible differences in storing of the floating point 

numbers in different programming languages or while using different compilers. 



13 



with the only elastic collision integrals, implemented in the calculations. There are also 
some additional arguments for the usage of these cross sections. From the phenomenological 
considerations one can take into account shortening of the mean free paths (or enlarging 
of the collision rates) due to creation of the resonances. In other words, there would have 
to be contributions from the inelastic 2 -H- 1 collision integrals next to the elastic collision 
integrals, and they are taken into account approximately by the contributions from the 
quasielastic cross sections. 

Resonances are not just intermediate particles, and they can collide with other particles. 
They make not negligibly small contribution to the thermodynamic functions and the vis- 
cosities. So that they are also included in the calculations as independent particles with their 
parameters and corresponding 2 -h- 2 collision integrals. They would have to have shorten- 
ing of their mean free paths from their decays and contributions from the inelastic 2 -H- 1 
collision integrals too. These contributions may be taken into account from the following 
collision rate considerations. A resonance's decay rate can be approximately replaced with 
just its total width. Then, given a resonance, one would have to redistribute its width (that 
is not changing the whole collision rate containing the contribution of the decay rate) in 
such a way that the cross section of the collision of this resonance with a resonance of the 
same species gets an addition 25 . Using an approximate expression for the collision rates 
(in the nonrelativistic approximation, applicable in this case) from Appendix |Dj one easily 
finds the addition r&/ (\/2nk{\vk\)) (where T k is the width) to the Aira^- Such cross sections 
seem to be the most physically preferable ones because they take into account more realistic 
mean free paths than in the previous case while not violating the conservation of the particle 
numbers too. 

The TCSs are used to take into account even larger shortening 26 of the mean free paths 
than in the case of the EQCSs. However, such cross sections introduce some inconsistency, 
implying that the conservation of the particle numbers is violated. As long as there may 
be contributions from some partial cross sections to the UrQMD ECSs or the EQCSs which 
were not taken into account (see below), the TCSs can be used as the upper bounds for the 
ECSs and the EQCSs. However, it's expected that these bounds are excessively high. If 
so, the TCSs (rather TCS2s) can be considered not only as the approximation taking into 
account real mean free paths but also as some measure of deviation from the approximation 
of only the elastic and the quasielastic collisions with the following arguments. If the TCSs 
were approximately equal to the ECSs or the EQCSs, or the numbers of particles with large 
inelastic cross sections were small, then one could expect small errors due to the negligibility 
of the inelastic collisions. 

Continuing the discussion of the details of the UrQMD cross sections, it should be men- 
tioned that the UrQMD TCSs are the most reliable ones. The sum of the partial cross 
sections is not always equal to the TCSs by their construction. If this is the case, then some 



This enhancement leads to the shortening of the mean free paths of the resonances of only this species, as 
needed. In the formal limit of this infinitely large enhancement other collision integrals can be neglected 
and the Boltzmann equation for this species decouples. Then from the solution of the Boltzmann equation 
for a single-component gas (see Sec. IVII A[) one concludes that the nonequilibrium perturbation to the 
distribution function of this species vanishes in this limit. Note that infinitely strong interactions also 
with particles of all other particle species would result in zero transport coefficients. 
Not the largest one. The effect of the enhancements of the resonances' TCSs is of 5% for the bulk viscosity 
and of 50% for the shear viscosity so that TCS2s are additionally considered. 
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partial cross sections are rescaled depending on their reliability 27 . 

The magnitudes of the partial cross sections, implemented in the UrQMD codes, are used 
to determine, what a partial cross section to choose in a given collision, using a random 
number generator. Among these partial cross sections there are the ECSs. Exactly these 
ECSs are used in the present calculations. However, if a partial cross section with a string 
excitation is chosen in a given collision, there is a probability to end up with the elastic 
collision if the yfs is too small. These contributions to the ECSs are not calculated and are 
not added to the ECSs. Also the string excitations can, possibly, end up with creation of a 
resonance. Contributions to the EQCSs from the string excitations are taken into account 
partially (see below). 

The ECSs, if not known from the experiment, are taken in the form of some extrapola- 
tions, discussed below, or the AQM is used. The normalization on the corresponding TCSs 
can change the ECSs notably. The meson meson (MM) ECSs are equal to 5 mb. The meson 
baryon (MB) ECSs are equal to the AQM rescaled experimental ir + p cross sections. But 
after the normalization they become equal to zero in the resonances dominated energy range 
(approximately below y/s = 1.7 GeV). The anti-baryon baryon (BB) ECSs are equal to the 
AQM rescaled experimental pp cross sections. Other ECSs are equal to the AQM ECSs. 

Before discussing the quasielastic cross sections first let's write for convenience the reso- 
nance dominated cross sections formula for a reaction 2 — > 1 — > any. Correcting a typo and 



rewriting it in a somewhat different form than in [30l . l3lj , one gets 



9r tt T 2 Rtot b R ^ij . Tj 

^ — , bn^ij = \{j t ,m i ,j j ,m j \\J R ,M R )\ — 

R ijinj rem v "/ ' ~ R,tot/ ^ *- R J ot 



(3.2) 

where is the partial energy-dependent width of the decay of the resonance R into 

particles of types i and j without specification of their isospin projection, T Ritot is the total 
energy-dependent width of the decay of the resonance R, gi is the spin degeneracy factor, 
b R -^ij is the energy- dependent branching ratio. The squared Clebsch-Gordan coefficients 
allow to specify the branching ratio b R -+ij for the pair of the particles with concrete isospin 
projections. The squared Clebsch-Gordan coefficients should be normalized in such a way 
that they give unity after summation over all isospin projections in a given multiplet. This 
formula represents contributions from all possible resonances through which the reaction can 
take place. Now it's easy to write down the cross sections for the quasielastic 2 — > 1 — > 2 
scatterings: 

ij ( /T\ _ ULJL T R,tot b R^ij /„ o\ 



R 



gi g jP l m (M R -^ + T% itot /4 



One more multiplier b R ->ij takes into account the fact that a resonance R decays only into 
the ij pair and represents the probability of this decay. 

The K~p TCS is not described by the formula (13. 2 h completely, and a partial cross 
section, attributed to the s-channel strings excitations, is added in the UrQMD model to fit 
the TCS to the experimental data. In the UrQMD model this s-channel strings cross section 



This information, including some other information about the cross sections, is stored in the array SigmaLn 
of the file blockres.f 
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is added also to other strange meson nonstrange baryon TCSs when annihilation is possible 
due to the quark content. From comparison with the experimental data for the K~p ECS 
*5l\ (actually it's believed to be the EQCS because smaller peaks from the resonances are 



reproduced there) it was found that the half of the s-channel strings cross section is enough 
to describe well this experimental K~p cross section. Then the half of the s-channel strings 
cross section is added to other strange meson nonstrange baryon EQCSs when annihilation 
is possible. These contributions from the strings excitations are the most low energetic ones. 
They are the only contributions from the strings excitations which are added. The next in 
the energy scale possible contributions to the EQCSs may be in the BB cross sections. In 
other pairs the string excitations appear approximately from yfs = 3 GeV . 

There is an important omission, found in the UrQMD codes (present also in the last 
version 3.3). The function fcgk returns incorrect (two times smaller) values of the squared 
Clebsch-Gordan coefficients for the resonances dominated cross sections in some cases. The 
first case is for the pairs of unflavored mesons from the same multiplet with the isospin 1 = 1. 
For example, the function fcgk returns 0.5 for the only possible isospin decomposition of the 
p + to the 7T + 7r° pair, because the states 7r + 7r° and 7r°7r + are counted as different ones. As 
a result, the peak from the p-resonance becomes two times smaller than e. g. in the p°- 
resonance isospin decomposition. The second less important case is for the pairs of unflavored 
mesons with the isospin 1 = 1 and anti-nucleons. The third even less important case is for 
the pair KK* and it's charge conjugate. 

Let's make some comments on the errors what the above-mentioned omissions cause in 
some quantities at zero chemical potentials, which in turn demonstrate sensitivity to different 
changes in the cross sections. The errors in the viscosities with the ECSs are less than 2%. 
The errors in the shear viscosity with the EQCSs (the TCSs) reach 57 — 63% (29 — 32%) at 
T = 0.07 GeV . Outside the temperature range 0.03 GeV <T< 0.14 GeV the errors reach 
11.6% (5.3%). The errors in the bulk viscosity with the EQCSs (the TCSs) reach 14.4-15.4% 
(10.6-11.4%) at T = 0.07 GeV. Outside the temperature range 0.03 GeV < T < 0.13 GeV 
the errors reach 4.8% (2.1%). The errors in the total number of collisions per unit time per 
unit volume (using the TCSs and including the decay rates) reach 10.2% (at T = 0.07 GeV). 
Outside the temperature range 0.04 GeV <T< 0.14 GeV the errors reach 5.1%. In view of 
the errors for the total number of collisions the kinetic freeze-out temperatures found in the 
UrQMD studies 49] should decrease, becoming closer to the experimentally extracted ones 
(see [58] and references therein). The chemical freeze-out temperature may change in a less 
extent. This is because both the inelastic and the quasielastic processes' cross sections (like 
of the quasielastic collision of the 7r + 7r° pair and of the reaction tt + ti — > K + K°) increase, 
so that the temperature at which the inelastic processes cease to be dominant may almost 
not change. 

It's observed that some of the UrQMD detailed balance cross sections (e. g. for the 
A + A° pair) are not symmetric under the particle interchange. This is because the function 
W3j, calculating the Wigner 3 — j symbols, doesn't return zero in some cases. Namely, the 
selection rule for the sum Ji + J 2 + J3 is not included. In principle, such omission could 
result in negative values of the essentially non- negative viscosities but, as long as only small 
fraction of cross sections is affected, this omission has caused only negligibly small errors in 
the viscosities. But e. g. the error in the A + A° TCS is approximately 25%. 

Also some fixes of the UrQMD cross sections are made. It's found that the 



K + p UrQMD ECS has large deviations from the experimental data [57J in the range 



1.6 GeV < \/s < 5 GeV (the UrQMD K + p cross section reaches 17 mb in the region 
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1.9 GeV < \fs < 3.1 GeV). To fit this cross section to the experimental data it is replaced 
with the AQM ECS in the range 3 GeV < \fs < 5 GeV and is interpolated smoothly with 
the sine function in the range 1.6 GeV < y/s < 3 GeV with the cross section being equal 
to 12.5 mb at \fs = 1.6 GeV. This replacement is also applied to other MB ECSs, when 
annihilation is not possible due to the quark content. 

The next fix is for the BB ECSs. It's found that the Ap UrQMD ECS has quite large 
deviations form the experimental data |57| too. To fit this cross section to the experimental 
data it is replaced with the AQM ECS in the range 4 GeV < y/s < 5 GeV and interpolated 
smoothly with the sine function in the range 2.2 GeV < \fs < 4 GeV with the cross section 
being equal to the AQM TCS at \fs = 2.2 GeV. This replacement is also applied to other 
BB ECSs. 

Some other found lacks result in negligible errors in the viscosities. However, errors in 
the corresponding mean free paths and possible other quantities may be not negligible ones. 
Two of such lacks can be mentioned. The first one is the following. The np and KN cross 
sections are fitted to the experimental data. And their charge conjugates are calculated 
using general formulas and so cause deviations up to 50% for \fs > 1.7 GeV . The second 
lack is the following. In some not large energy regions with \fs < 1.7 GeV the resonance 
dominated cross sections are equal to zero for some small numbers of pairs because there is 
no resonances which could be created by this pair. These regions are replaced by a constant 
continuously. 

Let's also comment on the deviations from the fixes described in the last four preceding 
paragraphs. The altogether deviations in the viscosities and the total number of collisions 
with the TCSs are less than 0.1%. The altogether deviations in the viscosities with the ECS 
or the EQCS are in the range 21 — 27%. The largest contribution is from the MB cross 
sections' fixes. At T < 0.12 GeV the deviations are less than 5% (the temperatures above 
T = 0.27 GeV are not studied). 

IV. CONDITIONS OF APPLICABILITY 

Before proceeding forth first the applicability of the Boltzmann equation and of the 
calculations of the transport coefficients should be clarified. 

Although the Boltzmann equations are valid for any perturbations of the distribution 
functions they should be slowly varying functions of the space-time coordinates to justify 
that they can be considered as functions of macroscopic quantities like the temperature, 
the chemical potentials or the flow velocity or, in other words, that one can apply ther- 
modynamics locally. Then one can make the expansion over the independent gradients of 
the thermodynamic functions and the flow velocity (the Chapman-Enskog method), which 
vanish in equilibrium. Smallness of these perturbations of the distribution functions in com- 
pare to their leading parts ensures the validity of this expansion and that the gradients 
are small 28 . Because these perturbations are inversely proportional to coupling constants 
one can say that they are proportional to some product of particles' mean free paths and 



The magnitudes of thermodynamic quantities can also be restricted by this condition or, conversely, not 
restricted even if transport coefficients diverge. See also Sec. IVII Al of this paper. The smallness of the 
shear and the bulk viscosity gradients can also be checked by the condition of smallness of the T^ 1 '^ 
(|5.16p in compare to the Of course, the next corrections should be small too. 
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the gradients. So that, in other words, the mean free paths should be much smaller than 
characteristic lengths, on which the macroscopic quantities change considerably 29 . 

As is discussed in Sec. HI1 the inelastic processes may need addition treatment in the 
calculations of the bulk viscosity. There is a need to specify reasonable conditions when the 
inelastic processes can be neglected. One could use the following reliable criterion, which 
takes into account both the particle number densities and the intensity of the interactions: 

ta rV(t) 

dt j d 3 x Yl R Pi<^ (4- 1 ) 

*! nG all channels 

where V(t) is the system's volume, J2 ne a n channels Rffit ls ^ ne number of reactions of particles 
of the fc'-th species 30 over all channels per unit time per unit volume (analog of fIDip ). t\ 
is chosen to satisfy the inequality, and t 2 is equal to the moment of time at which the 
divergence of the flow velocity is relaxed (if this time can be estimated reliably with remained 
inelastic processes) or to the moment of time at which the system becomes practically 
not interacting (after expansion) because of large cumulative mean free path in compare 
to the system's size. Though this criterion is likely to be too strict, and at some higher 
temperatures the approximation of conserved particle numbers should still work well. The 
main alternative criterion is based on comparison of collision rates of elastic and inelastic 
processes (as implemented in the j49|). Using this criterion and some other ones, the chemical 



freeze-out line 31 in the T — \i B plane can be built for the hadron gas, see e. g. [59|, |60 
At zero chemical potentials the chemical freeze-out temperature is approximately equal to 
T c h.f. = 160 — 170 MeV. The remaining question is how good is the approximation of only 
the elastic collisions at T < 160 MeV. From the hydrodynamical description of the elliptic 
flow at RHIC it's found that £/s < 0.05 near the chemical freeze-out [6l|. The constant 
value £/s = 0.04 provides a good description of the elliptic flow both at RHIC and LHC 



62J. It seems that the approximation of conserved particle numbers is not implemented in 
the bulk viscosity formula used in the [37[ ■ The bulk viscosity obtained from it is very close 
to the one of this paper. These results support the choice of the approximation of only the 
elastic collisions at T < 160 MeV and show that the deviations are likely no more than in 
2-3 times. Anyway, the numerical calculations by the Kubo formula through simulations 
of collisions are desirable along and around the chemical freeze-out line for more accurate 



calculations (though the procedure of collisions of particles introduce some errors itself [30 
which should be kept in mind). 

Errors due to the Maxwell-Boltzmann statistics, used instead of the Bose-Einstein or 
the Fermi-Dirac ones, were found to be small for the vanishing chemical potentials 32 . Ac- 



cording to calculations for the pion gas in [281 ]. the bulk viscosity becomes 25% larger at 



T = 120 MeV and 33% larger at T = 200 MeV for the vanishing chemical potential. Al- 
though the relative deviations of the thermodynamic quantities of the pion gas at the nonvan- 



It's clear that the mean free paths should be smaller than the system's size too. 

Primed indexes run over the particle species without regard to their spin states. This assignment is 
clarified more in Sec. IV Al 

This is an approximation. In fact this should be a range in which particles of different particle species 
have their own freeze-out points. 

It should be mentioned that if the particles of the &;-th particle species are bosons and if Hk{x^) > 
then there is a (local) Bose-Einstein condensation for them, which should be treated in a special way. 
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ishing chemical potential \i = 100 MeV are not more than 20% 33 the bulk viscosity becomes 
up to 2.5 times more. The shear viscosity becomes 15% less at T = 120 MeV and 25% less 
at T = 200 MeV for the vanishing chemical potential and 33% less at T = 120 MeV and 
67% less at T = 200 MeV for the fi = 100 MeV. The corrections to the bulk viscosity of 
the fermion gas, according to calculations of the bulk viscosity source term, not presented 
in this paper, are of the opposite sign and approximately of the same magnitude. So that 
for the hadron gas the error due to the used classical statistics can be even smaller than for 
the pion gas. 

The numerical calculations in Sec. I VI I of the viscosities with the total cross sections 
justify the choice of one constant cross section for all hadrons. It's approximately equal to 
20 mb, corresponding to the effective radius r = 0.4 fm (as given by the (13 .11) ). which is 
used in the estimations below. 

The condition of applicability of the ideal gas equation of state is controlled by the di- 
mensionless parameter vn which appears in the first correction from the binary collisions in 
the virial expansion and should be small. Here v = 167rr 3 /3 is the so called excluded volume 
parameter and 1/n is the mean volume per particle. One finds vn w 0.09 at T = 120 MeV, 
vn Rj 0.2 at T = 140 MeV and vn 1 at T = 180 MeV for the vanishing chemical poten- 
tials. Along the chemical freeze-out line (its parametrization can be found in [36|) the 
vn grows from 0.07 to 0.49 with the temperature. From comparison with lattice calcula- 



tions [63j one can find that the corrections to the ideal gas equation of state are small at 
T < 140 MeV. One could suspect that even small corrections to the thermodynamic quan- 
tities can result in large corrections for the bulk viscosity, though this seems to be not the 
case. The errors in the bulk viscosity from the scale- violating contributions of the hadrons' 
masses are less than the errors from the contributions to the trace of the energy-momentum 
tensor (for more details see Sec. IVT]) . 

One more important requirement, which one needs to justify the Boltzmann equation 
approach, is that the mean free time should be much larger than h/fl (the Q is the char- 
acteristic single-particle energy) [64[ or the de Broglie wavelength should be much smaller 
than the mean free path [25[ to distinguish independent acts of collisions and for particles 
to have well-defined on-shell energy and momentum. This condition gets badly satisfied for 
high temperatures or densities. The mean free path of the particle species k' is given by 
the formula (ID 1 7[) or the formula (ID12[) if the inelastic processes can be neglected. The 
wavelength can be written as \ k > ~ ^ , where the averaged modulus of the momentum of 

the k'-th species {\p k >\) is 

_ jd 3 p k ,\p k ,\fj°\p k ,) _ 2e-^T(3 + 3^ + 4) / 8m k ,T K 5/2 {z k ,) 

where z k > = my/T, K 2 (x) is the modified Bessel function of the second kind. As it follows 
from the H4.2[) the largest wavelength is for the lightest particles, the 7r-mesons. The elastic 
collision mean free paths are close to each other for all particle species. Hence, the smallest 
value of the ratio Xy/ly ls f° r the tt- mesons. Its value is close to the value of the vn and 
is exponentially suppressed for small temperatures too. At the temperature T = 140 MeV 



33 The relative deviations of the thermodynamic quantities grow with the temperature for some fixed value 
of the chemical potential and tend to some constant. 
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(180 MeV) and the vanishing chemical potentials this ratio is equal to 0.18 (0.7). Along the 
chemical freeze-out line it grows from 0.12 to 0.37 with the temperature. 

To go beyond these conditions one can use the Kubo (or Green-Kubo) formulas, for in- 
stance. In the [23| the Kubo formulas were used to perform perturbative calculations of the 
viscosities in the leading order. Basing on this result, an example of effective weakly cou- 
pled kinetic theory of quasiparticle excitations with thermal masses and thermal scattering 
amplitudes was presented in the j24[. There the function U(q) (appearing because of the 
temperature dependence of the mass) takes into account the next in the coupling constant 
correction to the energy-momentum tensor and the equation of state 34 . For further devel- 
opments see 
therein. 



25l 1431 . 144j |. For some other approaches see |65H67l| and [68J with references 



V. DETAILS OF CALCULATIONS 

A. The Boltzmann equation and its solution 

The calculations in this paper go close to the ones in [42[ though with some differences 
and generalizations. Let's start from some definitions. Multi-indices k,l,m,n will be used 
to denote particle species with certain spin states. Indexes k', V, m', nl will be used to denote 
particle species without regard to their spin states (and run from 1 to the number of the 
particle species N') and a, b to denote conserved quantum numbers 35 . Quantifiers V with 
respect to the indexes are omitted in the text where they may be needed which won't result 
in a confusion. Because nothing depends on spin variables one has for every sum over the 
multi-indexes 

J]- = (5.1) 

k k' 

where gy is the spin degeneracy factor. The following assignments will be used: 



n h = } n k i, n a = ) q ak n k , x k = — , x k > = — , x a = — , 
z — 4 — 4 n n n 

k k' k 

" = V± ?, - mk A* = Pk_ _ VV^IX s 

l^k — rp 1 H-a — rp ) Z k — rp ? 71 k ~ rp ) T k — rp > \°-^) 

where q ak denotes values of conserved quantum numbers of the a-th kind of the fc-th particle 
species. Everywhere the particle number densities are summed, the spin degeneracy factor g k i 
appears and then gets absorbed into the n k > or the x k > by the definition. All other quantities 
with primed and unprimed indexes don't differ, except for rates, the mean free times and 



34 In the hadron gas it's believed that the vacuum masses are large in compare to their thermal corrections 
for the most of the hadrons at temperatures T < 200 MeV or even higher ones. Then expanding over the 
thermal correction in the matrix elements, one would get even smaller corrections than the ones to the 
equation of state in coupling constants (because of coupling constants next to the matrix elements) in a 
perturbation theory, e. g. chiral perturbation theory. 

In systems with only the elastic collisions each particle species have their own " conserved quantum num- 
ber", equal to 1. 



20 



the mean free paths defined in Appendix [Dj the commented below, the coefficients A r k q,, 
Clf v and, of course, quantities, whose free indexes set the indexes of the particle number 
densities n^. Also the assignment J = J will be used for compactness somewhere. 

Pfz Pk 

The particle number flows are 36 

N ^S$kf* h ' (5 - 3) 

where the assignment fk{Pk) = fk is introduced. The energy- momentum tensor is 

^£/(J^PM/, (5-4) 

The local equilibrium distribution functions are 

fP = e^-^ T , (5.5) 

where is the chemical potential of the k-th particle species, T is the temperature and 
is the relativistic flow 4-velocity such that U li U >i = 1 (with a frequently used consequence 
U^dyU^ = 0). The local equilibrium is considered as perturbations of independent ther- 
modynamic variables and the flow velocity over a global equilibrium such that they can 
depend on the space-time coordinate . Additional chemical perturbations could also be 
considered, but they don't enter in the first order transport coefficients if they are small, 
as is discussed in Sec. [Til The chemical equilibrium implies that the particle number den- 
sities are equal to their global equilibrium values. The global equilibrium is called the 
time- independent stationary state with the maximal entropy 37 . The global equilibrium of 
an isolated system can be found by variation of the total nonequilibrium entropy functional 
'691 ] over the distribution functions with condition of the total energy and the total net 



charges conservation: 

W - £ / - -A) - £ / -£**,/ < 5 -«) 

where 0, X a are the Lagrange coefficients. Equating the first variation to zero, one easily 
gets the function (J5J| with = (1, 0, 0, 0), /3 = J and 



5^<M/ia, (5.7) 



where /i a = A a are the independent chemical potentials coupled to the conserved net charges. 

With f k = fl°\ substituted in the (ESJ) and the (Oil , one gets the leading contribution 
in the gradients expansion of the particle number flow and the energy-momentum tensor: 

= n k U", (5.8) 



The +,—,—,— metric signature is used throughout the paper. 

The kinetic equilibrium implies that the momentum distributions are the same as in the global equilibrium. 
Thus, a state of a system with both the pointwise (for the whole system) kinetic and the pointwise chemical 
equilibria is the global equilibrium. 
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where the projector 

A"" = 0"" - EW, (5.10) 
is introduced. The n k is the ideal gas particle number density, 

n k = UX 0) " = ^z\K 2 {z k )&\ (5.11) 
the e is the ideal gas energy density, 

e = W <»>- = £ / |^ P ;/(<» = £ n t e t , e k s m ,|||| - T, (5.12) 

and the P is the ideal gas pressure, 

p = -ir<»)-A^ = £ i / ^jfi'/f = X> r = (" 3 > 

Also the following assignments are used: 

e = - = x k e k , h k = e k + T, h = — — = x fc /i fc , 
n z — ' n L — ' 

k k 

e k K 3 (z k ) A _ e t h k K 3 (z k ) t h 

e k = — = z k - 1, e = -, h k = — = z k — h = -. 5.14 

T i^fc) T T K 2 (2fc) T 

Above h is the enthalpy per particle, e is the energy per particle and h k , e k are the enthalpy 
and the energy per particle of the fc-th particle species correspondingly, which are well 
defined in the ideal gas. 

In the relativistic hydrodynamics the flow velocity needs somewhat extended defini- 
tion. The most convenient condition which can be applied to the is the Landau-Lifshitz 
condition [4(| (Section 136). This condition states that in the local rest frame (where 
the flow velocity is zero though its gradient can have a nonzero value) each imaginary in- 
finitesimal cell of fluid should have zero momentum, and its energy density and the charge 
density should be related to other thermodynamic quantities through the equilibrium ther- 
modynamic relations (without a contribution of nonequilibrium dissipations). Its covariant 
mathematical formulation is 

{T» v - T (0) ^)[/ M = 0, {N£ - Nf^U^ = 0. (5.15) 

The next to leading correction over the gradients expansion to the T^ v can be written as an 
expansion over the 1-st order Lorentz covariant gradients, which are rotationally and space 
inversion invariant and satisfy the Landau-Lifshitz condition 38 (15.151) : 

T (i)„u = + £A pu V p U p = 7] ^A£A^ + ApA^ - j^A^A^ V P U T + £A^V P £/ P , 

(5.16) 



:i,s 



Also this form of T^ fJ,l/ respects the second law of thermodynamics [4Cf (Section 136). 
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where for any tensor a^ u the symmetrized traceless tensor assignment is introduced: 

( A ^ T ^ A ^ AtT " oT = A, upT a^, A^ pT A^ x = A^ A . (5.17) 



The equation ( 15 . 16f) is the definition of the shear t] and the bulk £ viscosity coefficients. The 
£A^V P £/ P term in the ( 15. 16ft can be considered as a nonequilibrium contribution to the 
pressure which enters in the (15. 9p . 

By means of the projector (I5.10P one can split the space-time derivative <9 M as 

= U^U u d u + = U„D + V„, (5.18) 

where D = U u d v , = A^d v . In the local rest frame (where = (1,0,0,0)) the D 
becomes the time derivative and the becomes the spacial derivative. Then the Boltzmann 
equations can be written in the form 

*W* = (Pk U » D + tfFM = Cf [A] + CrVtl (5.19) 

where C l k nel [fk\ represents the inelastic or number- changing collision integrals (it is omitted 
in calculations in this paper if the opposite is not stated explicitly) and C k [fk] is the elastic 
2f>2 collision integral. The collision integral Cj*[fk] has the form of the sum of positive 
gain terms and negative loss terms. Its explicit form is 39 (cf. (231. |25|) 

r ei in _ V- I f d3 Pu d3 P'k d3 P'u {rr _ t t \ 
k[Jkl Z-. lM 2 J (2vr) 3 2^ (2tt) 3 2^ (27r) 3 2< lJfcJl/ JkJll) 

x \M M \ 2 (27T) 4 6 4 (p' k +p' ll -p k - Pll ), (5.20) 

where jki = | if k and / denote the same particle species without regard to the spin states and 
■jki = 1 otherwise, \M.ki{p'kiPwPk-,Pii)\ 2 = l-M^I 2 is the square of the dimensionless elastic 
scattering amplitude averaged over the initial spin states and summed over the final ones. In- 
dex 1 designates that p k and p\ k are different variables. Introducing Wki = Wki(p' k , p'u',Pk,Pu) 
as 

W kl = l -^^5Vk + p'ii-Pk-Pii)i (5.21) 
one can rewrite the collision integral (I5.20p in the form as in [42j (Chap. I, Sec. 2) 

The Wki is related to the elastic differential cross section a k i as [42| (Chap. I, Sec. 2) 

W k i = sa k iS 4 (p' k + p' u - p k - pu), (5.23) 

where s = (p k + pu) 2 is the usual Mandelstam variable. The Wki has properties 
W k i{Pk,p'ii,Pk,Pii) = Wki(pk,Pu;p'k,Pu) = Wik(p' w p k ]PihPk) (due to time reversibility and 



39 The factor jf-l cancels double counting in integration over momentums of identical particles. The factor 
i comes from the relativistic normalization of the scattering amplitudes. 
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a freedom of relabelling of order numbers of particles taking part in reaction). And e.g. 
Wki(p' k ,p'u;Pk,Pn) 7^ Wki(p' w p' k ;pu,Pk) in the general case. The elastic collision integrals 
have important properties which one can easily prove [42| (Chap. II, Sec. 1): 

Swk, ci[h] = a - (5 ' 24) 



E/^J»1 = »' (5-25) 

Also the Cf[fk] vanishes if fk = fj®. 

The distribution functions fk solving the system of the Boltzmann equations approxi- 
mately are sought in the form 

fk = /f + f { k 1] = fl 0) + (5-26) 

where it's assumed that fk depend on the x^ entirely through the T, or their space- 

time derivatives. Also it is assumed that \(pk\ <C 1. After substitution of fk = in the 
(I5.19P the r. h. s. becomes zero and the 1. h. s. is zero only if the T, \x\. and U p don't depend 
on the x^ 1 (provided they don't depend on the momentum p k ). The 1-st order space-time 
derivatives of the T, fi k , U p in the 1. h. s. should be cancelled by the first nonvanishing 
contribution in the r. h. s. This means that the ipk should be proportional to the 1-st order 
space-time derivatives of the T, fi k , U^. The covariant time derivatives D can be expressed 
through the covariant spacial derivatives by means of approximate hydrodynamic equations, 
valid at the same order in the gradients expansion. Let's derive them. Integrating the (I5.19P 
over the ^p^o with the fk = in the 1. h. s. with the inelastic collision integrals retained 

and using the (15.241) and the (15.31) one would get (which can be justified using explicit form 
of the inelastic collision integrals) 

d^N^ = Dn k + n k V,U» = I k , (5.27) 

where I k is the sum of the inelastic collision integrals integrated over the momentum. It is 
responsible for the nonconservation of the total particle number of the k-th particle species 
and has the property q a kh = 0. If C l k nel [fk] = 0, then I k = which results in conservation 
of the total particle numbers of each particle species. Multiplying the (I5.27P on the q ak and 
summing over k one gets the continuity equations for the net charge flows: 

d„N^ = Dn a + n a VJJ» = 0. (5.28) 
Then integrating the (I5.19P over the p k ^y^p with the fk = in the 1. h. s. one gets 

d pT (0)p» = d p (eU p U v - PA pv ) = 0. (5.29) 

There is zero in the r. h. s. even if the inelastic collision integrals are retained because 
they respect energy conservation too. Note that the Boltzmann equations (I5.19P (without 
any thermal corrections) permit a self-consistent description only if the energy-momentum 
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tensor and the net charge flows of the ideal gas are used. After the convolution of the (I5.29p 
with the one gets the Euler's equation: 

DW = —^—V^P = -!-V M P. (5.30) 
e + P hn v ' 

After the convolution of the (15.291) with the U u one gets equation for the energy density: 

De = -(e + P)V li U ,l = -hnV li U' i . (5.31) 

To proceed farther one needs to expand the 1. h. s. of the Boltzmann equations ( 15. 19)) 
over the gradients of thermodynamic variables and the flow velocity. Let's choose the \x a 
and the T as the independent thermodynamic variables. Then for the Df% one can write 
the expansion 

f>f(°) r)f(°) /}f(°) 

D ^=E^t D ^ + %r DT + %rDU>-. (5.32) 

Writing the expansion for the Dn a and the De one gets from the (I5.28P and the (15.311) : 

Dn a = J2 p^Df* + ^f DT = -n a V^, (5.33) 

b 

De = ^DT + J2 J~ D ^ = -hnVJJ". (5.34) 

The solution to the system of equations (15.331) . (15.34)) can be found easily: 

DT = -RTV^U", (5.35) 

Dfi a = Tj2 Kb (R B b - x h )VJJ», (5.36) 
b 

where 

R = ~ h -Z^»E^t± , (5 .37) 

and 

Oti ti ~ dfi ti de de 

df b =T Aab > W = f Ba > &T= nC ^ Wa^ nEa - (5 - 38) 

Above it is assumed that the matrix A ab is not degenerate 40 , which is related to the self- 
consistency of the statistical description of the system. Using the ideal gas formulas (15.1 ip 



40 One can prove that the N" x N" matrix A a b in (|5.39[) is not degenerate if there are N" linearly independent 
conserved charges. Then one can prove that the denominator in the (|5.40[) is not zero. 
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and (I5.12p one gets 

A ab = qakqbkXk, E a = ^2 Qak x kh, B a = E a - Aj&/*6, (5.39) 

k k b 

C M = ^2 Xk ^hk + 4~ h^k) = Xk(3h k + Z 2 k ) - 22 E aP-a = C {fl} - #aAa, 
fc fc a a 

and simplified expressions for the R and the £)/t a 

fl= J^iE^l , (5.40) 

~ Ea,fe Art 

L>/2 a = Yl A ab( RE b - ^) V M ^. (5.41) 
b 

For the special case of the vanishing chemical potentials, /i a — > 0, (for a chargeless system the 
result is the same) the quantities n a , x a , B a , E a tend to zero because the contributions from 
particles and anti-particles cancel each other and the chargeless particles don't contribute. 
Then from the (15.351) and the (15.361) one finds 

DT\, a =o = -J-W, (5.42) 

£>//aL=o = 0. (5.43) 

This means that for the vanishing chemical potentials one can simply exclude them from the 
distribution functions (if one does not study diffusion or thermal conductivity). In systems 
with only the elastic collisions each particle has its own charge so that one takes q a k = 5 ak 
and gets 

Aki = SkiXk, B k = x k (e k — ftk), E k = e k x k , R = — , 

c v 

C M - J2 E ^ab B b = J2%k(-ht + 5h k + z\ - 1) = J2x k c V:k = c v . (5.44) 

a, b k k 

Then the equation for the DT (I5.35P remains the same with a new R from the (I5.44p . and 
the equations (I5.36P become 

Dfi k = (J-{e k - fi k ) ~ T^j V^. (5.45) 

Note that in systems with only the elastic collisions the D[i k does not tend to zero for 
the vanishing chemical potentials so that the fi k could not be omitted in the distribution 
functions in this case. Because the heat conductivity and diffusion are not considered in 
this paper their nonequilibrium gradients are taken equal to zero, V U P = V V T = V u fi a = 0. 
Using the (15.351) . (15.361) and (15.301) the 1. h. s. of the (15.191) can be transformed as 

WJ„P + ^V M )/i 0) = -T/f + Tfl 0) Q k V P U", (5.46) 
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where 

Qk = r 2 k (l-R)+T k J2 QakA- b \RE b - x b ) - l -z\. (5.47) 

^ ' a,b 



Using the (15.171) one can notice that the useful equality vr^7r^ V fJJ v = Tc k 7r k V ^U v holds. In 
systems with only the elastic collisions the Q k simplifies in agreement with |42j (Chap. V, 
Sec. 1): 

where the assignments 7 from |42] is used. It can be expressed through the c v , defined in 
the (15.441) . as 7 = — + 1. Introducing symmetric round brackets 

( F - G )' s 4^k)^// fe)Gfe)e " T '- (5 - 49) 

and assignments 

o 

o 

a r k = (Q k ,Tl), 7fc = i^KK^k^ku), a r k = (l,T r k ) k , (5.50) 

and using explicit expressions of the a k from Appendix [A] one finds for the a k and the a\ in 
systems with elastic and inelastic collisions 



«2 



1 + Yl VakA- b \RE b - x b ) - e k R } (5.51) 



a,b 



ul = h k + Yl ^qakKl (R E b - x b ) - (3h k + z 2 k )R. (5.52) 

a.b 

Then using the (I5.5ip and the (I5.52p one can show that 

^2q a kX k a° k = 0, (5.53) 

k 

J2x k al = 0. (5.54) 



Because the gradients V M f/ M and V \JJ V are independent the (I5.53f) and the (15.54ft are direct 
consequences of the local net charge (I5.28P and the energy-momentum (I5.29P conservations. 

o o 

Quantities (1, vr^7r^) and (p k , Tt k Tt k ) vanish automatically because of the special tensorial 

o 

structure 41 of the it k 'K k . 



41 Direct computation gives (l.TrJftQ* oc {dVU" + C 2 A CT ")A^ p = 0, (p^.Tr^)* oc {dU x U a UP 
C 2 U x A' y P + C 3 U a A^A^ = 0. 
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The next step is to transform the r. h. s. of the Boltzmann equations (I5.19p . After the 
substitution of the (I5.26P in the r. h. s. of the (I5.19P the collision integral becomes linear 
and one gets 

^[/^-/1 0) E £ «N' ( 5 - 55 ) 
i 

where 

= 7^3 / fS\v>k + Vii ~ V'k ~ riiWu- (5-56) 

\ Z7T ) J P u,p' k ,p'u 

The unknown functions tp k are sought in the form 

= [~A k [p k )V ^ + C k (p k )rfffiVJJ^ , (5.57) 

where cr(T) is some formal averaged cross section, used to come to dimensionless quantities. 

o 

Then using the (I5.46P and the (I5.55p . and the fact that the gradients V^U^ 1 and V M {7„ are 
independent, the Boltzmann equations can be written as independent integral equations: 

4 = ^>^[M (5.58) 



i 

where the dimensionless collision integrals are introduced: 

I&bCk) = -7^£&\Xk). (5-60) 
niTcr(T) 

In the case of present inelastic processes the 1. h. s. of the (I5.58P is set by the source term 
( I5.47P and the r. h. s. contains the linear inelastic collision integrals. After introduction 
of inelastic processes the source terms in the ( I5.58P become much larger as demonstrated in 
Sec. IVH Al Using the equations (15.361) and (I5.35j) and the ideal gas formulas (15.391) one can 
check that in the zero masses limit the source terms Qk ( I5.47P tend to zero and Dfi a = that 
is the fi a don't scale and the distribution functions become scale invariant. The source term 
of the shear viscosity in the (15.59j) doesn't depend on the presence of inelastic processes in 
the system and originates from the free propagation term in the Boltzmann equation. 

B. The transport coefficients and their properties 

After substitution of the with the (p k ( I5.57P into the (15 .4p and comparison with the 
(I5.16P one finds the formula for the bulk viscosity 

1 T . 

£ = ~^~rr\ Xki^KukKvk, A k )k, (5.61) 
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and for the shear viscosity 



1 T x - - o 

^ = TT^^a MKK' C k^k^ku)k, (5-62) 

where the relation A^^A^A^ = 5 is used. 

In kinetics the conditions that the nonequilibrium perturbations of the distribution func- 
tions does not contribute to the net charge and the energy-momentum densities are used as 
a convenient choice and are called matching conditions. They reproduce the Landau-Lifshitz 
condition (I5.15p . The matching conditions for the net charge densities can be written as 

S^/(2TM^ V = ' (5 ' 63) 
and for the energy-momentum density can be written as 

For the special tensorial functions C k iT k ^7r ku i n the (I5.57P they are satisfied automatically 
and for the scalar functions A k they can be rewritten in the form (the 3-vector part of the 
(I5.64p is automatically satisfied) 

^qakXkin, A k ) k = 0, ^x fc (r^, A k ) k = 0. (5.65) 

k k 

The conditions (I5.63f) and (15.641) exclude the nonphysical solutions A z k - m - = J2 a G a Qak + Cr k 
(which cannot be solutions in inhomogeneous systems and are produced just due to shifts 
in the T, /z a ) of the linearized Boltzmann equations for which the collision integrals vanish 
(Afc™" are zero modes). From the formula (15.611) one can see that in the framework of the 
Boltzmann equation these conditions also eliminate ambiguity in the £ due to freedom of 
addition of the A k m ' to the solution A k of the (15.581) . With help of these matching conditions 
one can show explicitly essential positiveness of the £. Namely, using the matching conditions 
(I5.65p . the equation (15.581) and the identity A^n^ tk ir U;k = z \~ r f , the bulk viscosity (I5.6ip 
can be rewritten as 

£=^E^>^ = ^E Xfc (^ =^jr ) i{A},{A}}, (5.66) 

where the square brackets are introduced for sets of equal lengths {F} = (Fi, F k , ...), 
{G} = (Gi, G k , ...): 

{G}] = -^=r E 7^6 / fi 0) f$\Fk + F U - F' k - F' u )G k W M . (5.67) 

^ > k,l ^ ' " 'Pk'Pll'P 'k,P 'll 

Using the time reversibility property of the W k i one can show that the equality 

(F k + F u - F' k - F' u )G k = ±(F k + F ll - F' k - F' u )(G k + G u - G' k - G'u), (5.68) 
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holds under the integration and the summation in the (15.671) . Then one gets the direct 
consequence 

[{F}, {G}] = [{G}, {F}], [{F}, {F}] > 0. (5.69) 

This proves the essential positiveness of the £. Similarly using the (I5.59p . the shear viscosity 
can be rewritten in essentially positive form 



1 T 

V 



J2 Xk {^kK^k^k^k^j = ^^7j^^2 X k \^Xl L tl[Ck^Z],C k Tt k ^ kv 



10a(T) 

v ' k x /fc v ' k \ I / k 

1 ' [{Cwk^},{CW^Z}}. (5.70) 



10 <r(T) 



The considered variational method allows to find an approximate solution of the integral 
equations (I5.58P and (I5.59P in the form of a linear combination of test-functions. The 
coefficients next to the test-functions are found from the condition to deliver extremum to 



some functional. One could take this functional in the form of some special norm, as in |42 



Or one can take somewhat different functional, like in 70 , which is more convenient, and 



get the same result. This generalized functional can be written in the form 

F[ X ) = J2 x k(Sr,Xk»...»)k - \[{X"- U }, {X,...u}], (5.71) 



where S£" u = Q k and Xk^...u = A k for the bulk viscosity and S^'" v = Tr k n k , x k " V = CkK k 7r k 
for the shear viscosity. Equating to zero the first variation of the (15.711) over the Xkn...v one 
gets 

Y,Xk{S^\8 Xkl ,..,)k - [{X^},^...,}} = 0. (5.72) 

k 

Because the variations Sxkfi...u are arbitrary and independent the generalized integral equa- 
tions follows then: 

l 

The second variation of the (I5.7ip is 

S 2 F[ X ] = -[{5x' 1 - u },{Sx,...,}]<0, (5.74) 

which means that the solution of the integral equations (15.581) and (15.591) is reduced to the 
variational problem of finding of the maximum of the functional (15. 71 p . Using the (15.731) . 
the maximal value of the (I5.7ip can be written as 



:[x) = t;[{x' 1 - u }Ax,...Mx=x^- (5-75) 



Then using the (15.661) and the (I5.70p one can write the bulk and the shear viscosities through 
the maximal value of the F [x] 



T 

S = 2 T=rF„ 

a{T) 



(5.76) 
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This means that the precise solution of the (15. 73j) delivers the maximal values for the trans- 
port coefficients. 

The approximate solution of the system of the integral equations (I5.58P and ( I5.59P are 
sought in the form 

^ = E^ ( 5 - 78 ) 

r=0 



n 2 



Cu = ^Clrl (5.79) 

r=0 



where ri\ and n 2 set the number of the used test-functions. Test-functions, used in [70j . 
would cause less significant digit cancellation in numerical calculations but there is a need 
to reduce the dimension of the 12-dimensional integrals from these test-functions as more as 
possible to perform the calculations in a reasonable time. The test-functions in the form of 
just powers of the Tk seem to be the most convenient for this purpose. Questions concerning 
the uniqueness and the existence of the solution and the convergence of the approximate 
solution to the precise one are covered in (42] (Chap. IX, Sec. 1-2). As long as particles of 
the same particle species but with different spin states are undistinguishable their functions 
ifk (I5.57P are equal, and the variational problem is reduced to the variation of the coefficients 
A r k , and C£,, and the bulk (I5.66P and the shear (15.701) viscosities can be rewritten as 

rp N' m 



a(T) 

v ' k'=l r=0 



N' n 2 
^ ' k'=l r=0 

After the substitution of the approximate functions Ay (I5.78P and (I5.79P into the ( I5.7ip 

and equating the first variation of the functional to zero one gets the following matrix 
equations (with the multi-indexes (/', s) and (k',r)) for the bulk and the shear viscosities 
correspondingly 42 

N' m 

^<' = EE A ^, (5-82) 

l'=l s=0 



One can first derive the same equations for the Ak and Ck, treating them as different functions for all 
k, with the coefficients A r ^ and C r k \ having the same form as the A r ^i V and C£? ; , . Then after summation 
of the equations over spin states of identical particles and taking Ak — Ak> , Ck = Ck' one reproduces the 
system of equations for the Ak> and C'k> . 
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N> n 2 

^7& = EE C S'C?. ( 5 - 83 ) 

i'=l s=0 

where the introduced coefficients A T k , v and C k f v are 

iV' 

Ay v = X k >Xli[T r , T*] k >V + Sk'l'Xk' 2^ X m'V T i T %'m>, (5.84) 

m'=l 



AT' 



CfcT = Xk'Xi'l^nVn", TiW^7T^)k'i' + 2^ x m \T r tx^tx v , r s T^WZ\ klml . (5.85) 



m'=l 



They are expressed through the collision brackets 

[F, Gi] kl = —r- 2 . ^ . . -- f e-^ 1 (F k - F' k )G u W kl . (5.86) 



TeW4zfK 2 (z k )K 2 ( Zl )<T{T) 

The collision brackets [F, G] k i are obtained from the last formula by the replacement of the 
Gu on the G k . Due to the time reversibility property of the W k i one can replace the Gu on 
the \{Gn — G'ii) in the (15.861) . Then one can see that 

[r r : T s ] kl >Q. (5.87) 

Also it's easy to notice the following symmetries 

[F, Gi] kl = [G, Fi] lk , [F, G] k i = [G, F] kl . (5.88) 

They result in the following symmetric properties A r k s , v = A s J k ,, C T k f v = Cf^,. Also the micro- 
scopical particle number and energy conservation laws imply for the A^r. 

A% v = 0, (5.89) 



N' 

E A »v = °- ( 5 - 9 °) 

k'=l 

The (15.891) together with the a® = (1A13|) means that the equations with r = in the (I5.82p 
are excluded. From the (15.901) and (I5.54j) it follows that each one equation with r = 1 in the 
(I5.82p can be expressed through the sum of the other ones, reducing the rank of the matrix 
on 1. To solve the matrix equation (I5.82j) one eliminates one equation, for example with 
k' — 1, r = 1. One of coefficients of A], is independent; for example, let it be A\,. Using the 
(I5.90p . the matrix equation (15.821) can be rewritten as 

N' N' m 

x k ,al, = ^Alr k ,(Al - A\,) + ^^«. (5.91) 

V=2 l'=l s=2 
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Then, using the (1A13j) and the (I5.54p . the bulk viscosity (15 .80j) becomes 



T N' T N> ni 

^ > k'=2 ^ ' k'=l r=2 

Then the coefficient A\, can be eliminated by shift of other A\, and be implicitly used 
to satisfy one energy conservation matching condition. The particle number conservation 
matching conditions are implicitly satisfied by means of the coefficients A° k ,. The first term 
in the (I5.92p is present only in mixtures. That's why it is small in gases with close to 
each other masses of particles of different species (like in the pion gas). In gases with very 
different masses (like in the hadron gas) contribution of the first term in the (I5.92p can 
become dominant. 

Analytical expressions for some lowest orders collision brackets, which enter in the matrix 
equations ( I5.83P and (I5.9ip . can be found in Appendix Higher orders are not presented 
because of their bulky form. 



VI. THE NUMERICAL CALCULATIONS 

The numerical calculations for the hadron gas involve roughly 2 ^ N ^ the 12- dimensional 
integrals, where N' is the number of particle species and n is the number of the used test- 
functions (called the order of the calculations). The 12-dimensional integrals (they are 
the collision brackets [F, Gi]jy and [F, G]ki) can be reduced to 1-dimensional integrals. For 
constant cross sections they are expressed through special functions, and for other ones 
numerical methods are used. The details of calculations are described in Appendix This 
allows to perform the calculations with a good precision in a reasonable time. Because the 



is 



analytical expressions for the collision brackets are bulky the Mathematica package [71 
used for symbolical and some numerical manipulations. 

The calculations of the viscosities are quite reliable at T < 120 — 140 MeV (throughout 
the paper the chemical potentials are equal to zero if else is not stated), as is discussed in 
Sec. IIV1 The numerical calculations are done also for temperatures up to T = 270 MeV 
for the future comparisons and to show the position of the maximum of the bulk viscosity, 
when it is present only due to the hadrons' masses. Introduction of the inelastic processes 
should increase the bulk viscosity, though in the approximation when its nonzero value is 
maintained only by the hadrons' masses the maximum may shift not considerably. Taking 
into account the non-ideal gas equation of state, the maximum may shift to some extent too 
and become sharper, as can be seen from the speed of sound of the [HI, or a new smaller 
maximum can appear. 

The UrQMD (version 1.3) particle list is used, which doesn't contain charmed and bot- 
tomed particles and consists of 322 particle species including anti-particles. Some thermo- 
dynamical quantities of the ideal hadron gas with this mass spectrum are shown in fig. [2] 
and fig. [3j The e and P are given by the (I5.12p . (15. lip . (I5.13p . and the s is given by 
(1B7P . The quantity R (I5.40p , appearing in the bulk viscosity source term (I5.47P and tending 
to 1/3 if the particles' masses are tended to zero, is equal to the squared speed of sound 
Rch.-n. = v% = ^ in the case of the charge- neutral (or chargeless) system (implying equal 
to zero and not developing chemical potentials). The R is given by the R e i. c . (I5.44p in the 
case of only the elastic collisions. The R c h.- n . is quite close to the squared speed of sound 
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of the ideal gas in the [191 ] . For the intermediate case when there are conserved and not 
conserved particle numbers it's observed that the R is above 1/3. 

The new particle list with charmed and bottomed particles 43 (cut on 3 GeV, which 
results in negligible errors 0.01% or less) from the THERMUS package [73[ was used in the 



calculations at zero chemical potentials in the [26] . The errors in the trace of the energy- 
momentum tensor because of neglected rhe charmed and bottomed particles grow with the 
temperature and are equal to 13% (21%) at T = 140 MeV (270 MeV). The errors in the 
R are less than 3.2%. The errors in the shear viscosity (calculated with one constant cross 
section for all hadrons) are less than 0.6% and in the bulk viscosity are less than 7.5%. An 
additional study of the mass spectrum dependence of the viscosities can be found in the 
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FIG. 2: Some scaled thermodynamic quantities (the entropy density s, the energy density e, the 
pressure P) and the trace of the energy-momentum tensor of the ideal gas with the hadron gas 
mass spectrum including 322 particle species. 



The results for the shear and the bulk viscosities are shown in fig. H] and fig. [5] corre- 
spondingly. They are calculated with the different cross sections (with all corrections and 
fixes): ECQs, EQCSs, EQCS2s, TCSs, TCS2s. They are described in Sec. HH The bulk 
viscosity with the TCS2s is not shown because it's smaller only on 5% or less than the 
one with the TCSs. Up to 5 (3) test-functions are used in the calculations of the bulk 
(shear) viscosity. The maximal errors are 11% and are less than 4.2% outside the range 
40 MeV < T < 90 MeV. The best convergence is for the case of the ECSs (the errors are 
less than 2%). For quantitative results the recommended cross sections are the EQCS2s, as 
is commented in Sec. IIHI For qualitative analysis the TCSs and the TCS2s are more suit- 
able. Also it's shown that the approximation of one constant cross section for all hadrons 
is a good one. For the shear viscosity with the EQCSs or the EQCS2s it's somewhat worse. 
This can be explained in the following way. There are descending and growing cross sec- 
tions, and these opposite dependencies compensate approximately each other, so that at 
some temperatures the cross sections manifest themselves approximately as a constant one. 
Some EQCSs and EQCS2s have quite steeply descending tails, which explains relatively fast 



These are the particles which are more or less reliably detected [7! 
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FIG. 3: The quantity R (|5.40p . appearing in the bulk viscosity source term (|5.47p . ate zero chemical 
potentials. In the case of the charge- neutral (or chargeless) system the R coincides with the squared 
speed of sound v% = In the case of only the elastic collisions it is given by the (I5.44|) . It 
is calculated for the ideal hadron gas in the first and the second cases (solid and dashed lines 
correspondingly) and the ideal pion gas (doted and dash-dotted lines correspondingly). The value 
of 1/3 is shown for convenience (dash-dot-dotted line). 

rises in the shear viscosities, because of which one constant cross section doesn't provides 
worse approximation. An explanation why the bulk viscosity is approximated well in the 
same case with one constant cross sections at some temperatures would be somewhat more 
complicated. 

Also the calculations without resonances (the particles with the width larger or equal to 
0.2 MeV) are done with the EQCSs to find out the magnitude of their contributions. After 
the exclusion 26 particle species remain. The bulk viscosity decreases not more than in 2.8 
times, and the shear viscosity decreases not more than in 1.6 times (using the TCSs these 
factors are somewhat smaller). 

Note that at T pa 160 MeV the viscosities calculated with the TCS2s are approximately 
2 times smaller than the viscosities calculated with the EQCS2s respectively. This reflects 
the fact that the contribution to the total number of collisions from the inelastic processes is 
approximately the same as from the elastic plus the quasielastic processes at the freeze-out 
temperature. 

The maximum of the bulk viscosity is, of course, sensitive to the energy dependence of 
the cross sections, as can be seen from fig. E. g. if the BB and some MB EQCSs were 
not fixed, as is described in Sec. IHIt the maximum would be present at T pa 190 MeV. 
After the fixes the BB and some MB EQCSs and EQCS2s have steeply descending tails, 
and the maximum shifts towards much higher temperatures. Though for the qualitative 
analysis the TCSs and the TCS2s are more suitable because at T > 160 MeV the inelastic 
processes make not small contributions. With these cross sections the bulk viscosity has the 
maximum at T w 190 MeV. With one constant cross section the maximum is present at 
T pa 200 MeV with [26] or without charmed and bottomed particles. 

In several papers the viscosities of the hadron gas were studied by the ones of the pion 
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FIG. 4: The shear viscosity of the hadron gas, calculated with the ECSs, EQCSs, TCSs, TCS2s 
(dashed line) and the EQCS2s (solid line). One constant cross section for all hadrons is used in 
the approximating calculations. Different values are chosen: 5 mb, 8 mb, 20 mb (dotted line). See 
text for more details. 




FIG. 5: The bulk viscosity of the hadron gas, calculated with the ECSs, EQCSs, TCSs (dashed 
line) and the EQCS2s (solid line). One constant cross section for all hadrons is used in the 
approximating calculations. Different values are chosen: 5 mb, 8 mb, 20 mb (dotted line). See text 
for more details. 



gas. This approximations turn out to be bad while calculating the bulk viscosity 44 . In 
fig. [6] the ratio of the bulk viscosity of the hadron gas to the one of the pion gas, using 
different cross sections, is shown. With all the used cross sections the deviations are quite 
large. Purely pion gas implies the ECSs. The bulk viscosity of the hadron gas with the 
EQCS2s is divided on the one of the pion gas with the EQCS. This ratio reaches 122 at 



44 This discrepancy could be noticed earlier the results of the 37] and the 19(. Though they required 
confirmations or justifications. 
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T = 270 MeV. Considering a closer to the pion gas approximation, excluding the resonances 
in the hadron gas, this ratio becomes 2.3 — 2.6 times smaller at T = 120 — 140 MeV. The 
same factor is equal to 2.0 — 2.1 for the used TCS2s and the TCSs for the hadron gas and 
the pion gas correspondingly. Though the resonances should not be excluded. At the same 
time the corresponding ratios of the particle number densities and the energy densities at 
T = 140 MeV are approximately equal to 2 and 3 respectively. Note that the ratios of 
the viscosities are less sensitive to the corrections discussed in Sec. [IV] than the viscosities 
themselves. The ratio of the shear viscosity of the hadron gas to the one of the pion gas 
is not larger than 1.6, as can be seen in fig. [71 It can even be somewhat smaller than 1 
if the TCS2s and the TCS2 are used. This seems to be because the shear viscosity is not 
much sensitive to the mass spectrum as the bulk viscosity, and the contributions from the 
hadrons other than pions at high temperatures (where pion numbers don't dominate) come 
with somewhat larger (on average) cross sections. 

Some simplified explanations of the enlarged bulk viscosity in the hadron gas (to some 
extent) and the position of its maximum can be made. The bulk viscosity is sensitive 
to particle's masses m ~ T, and the hadron gas mass spectrum provides such masses at 
different temperatures. For a fixed m/T and approximately constant cross section the bulk 
viscosity grows with the temperature, as can be seen from the formula (17.41) . using it as an 
estimate. Particles with very large masses have relatively small number densities because 
of the exponential suppression e~ m l T so that the maximum is set by particles with not the 
largest masses at T 190 MeV, and a relatively slow further descending follows at higher 
temperatures. 




FIG. 6: The ratio of the bulk viscosity of the hadron gas to the one of the pion gas. The different 
cross sections are used: TCS2s/TCSs (dashed line), EQCS2s/EQCSs (solid line), ECSs/ECSs 
(dotted line). See text for more details. 



The ratio of the shear viscosity to the entropy density rj/s and the ratio of the bulk 
viscosity to the entropy density £/s in the hadron gas are shown in fig. El The EQCS2s 
are used. As long as the maximum of the bulk viscosity is not sharp, the ratio £/s doesn't 
have a maximum and is a descending function of the temperature. The entropy density is 
calculated by the formula (1B7|) using the ideal gas formulas in the (15. lip and the (15.141) . 
The ratio of the bulk viscosity to the shear viscosity with the EQCS2s is shown in fig. [91 



37 



1.6 




.4 1 ■ ' ■ ' ■ ' ■ "— 

0.05 0.10 0.15 0.20 0.25 

T (GeV) 



FIG. 7: The ratio of the shear viscosity of the hadron gas to the one of the pion gas. The different 
cross sections are used: TCS2s/TCSs (dashed line), EQCS2s/EQCSs (solid line), ECSs/ECSs 
(dotted line). See text for more details. 
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FIG. 8: The ratio of the shear viscosity to the entropy density and the ratio of the bulk viscosity 
times 100 to the entropy density in the hadron gas. The EQCS2s are used. 



The dependencies from the temperature of the r]/s and the £/s calculated along the 
chemical freeze-out line are found too and are depicted in fig. [TOj The EQCS2s are used. As 
was discussed in Sec. HV] the calculations with large chemical potentials may contain large 
deviations, especially in the bulk viscosity, however, in the hadron gas the contributions 
from the bosons and the fermions may cancel substantially. The calculations along the 
chemical freeze-out line could have not small deviations for the bulk viscosity because of 
the inelastic processes. At the considered collision energies the strange particle numbers 
are not described well by the equilibrium statistical calculations. It's expected that this is 
because they don't reach the chemical equilibrium before the chemical freeze-out takes place. 



After the introduction of the strange saturation factors 7^ |74[ the experimental data gets 
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FIG. 9: The ratio of the bulk viscosity times 100 to the shear viscosity in the hadron gas. The 
EQCS2s are used. 



described better. As long as the considered chemical perturbations are not quite accurate 
and are not small they are used in a phenomenological way, being inserted into all particle 
number densities 45 . Because of all this the calculations along the chemical freeze-out line are 
less reliable than the ones at zero chemical potentials. All variables' values of the chemical 
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FIG. 10: The ratio of the bulk viscosity times 100 to the entropy density and the ratio of the shear 
viscosity to the entropy density as functions of the freeze-out values of the temperature T. All 
other freeze-out line variables' values can be found in the [36]. The EQCS2s are used. 



The noncquilibrium chemical potential-like perturbations of the form Tn St k In 7 S (n Si k is the number of 
strange quarks in hadrons of the fc-th species) [7J], obtained from statistical description and reflecting 
suppression of the strange particle numbers, obviously violate conservation laws to some extent. However, 
they are used because of the simplicity in phenomenological estimating calculations. There are also more 
elaborated calculations of the chemical perturbations, see e. g. [75j . Also 



see 



581 ] for some discussions. 
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freeze-out line, including the strangeness saturation factor 7^, are conveniently presented in 
the [36]. The convergence of the calculations (with all the cross section types) is good with 



the errors less than 4% for both the viscosities. Also it was checked how the results change 
if the chemical equilibrium is assumed. The entropy density increases up to 36%. The shear 
viscosity increases no more than on 13%. The bulk viscosity decreases no more than on 
5.1% (though with the ECSs the decrease would be on 44%). 

VII. ANALYTICAL RESULTS 
A. The single-component gas 

In the single-component gas, using one test-function, the matrix equations can be easily 
solved, and the shear (15.811) and the bulk (I5.92p viscosities become (indexes "1" of the 
particle species are omitted) 

1 T (7 ) 2 

V = 10 ofT) ( 7 - X ) 
T (a 2 ) 2 

In this approximation the explicit closed-form (expressed through special and elementary 
functions) relativistic formulas for the bulk and the shear viscosities were obtained in the 
*39| . There the parameter a = 2r. In [42| (Chap. XI, Sec. 1) they are written through the 



parameter a = 2r 2 . 46 The results are 

15 T z 2 K 2 (z)h 2 



64tt r 2 (5z 2 + 2)K 2 {2z) + (3z 3 + A9z)K 3 (2z) 

= 1 T z 2 K 2 2 (z)[(5-3 7 )/*-3 7 ] 2 
5 64vrr 2 2K 2 (2z) + zK 3 (2z) ' 



(7.3) 
(7.4) 



where 7 = j- + 1 = J^t^^^ - Though the correct result for the shear viscosity is 

= l^T z 2 K 2 {z)h 2 

77 64nr 2 (15z 2 + 2)K 2 (2z) + (3z 3 + 49z)K 3 (2z)' { ' } 

This result is in agreement with the result in [27|, 76]. To get the ( I7.5P and the f)7.4p the 
collision brackets in the C 00 (I5.85P and the A 22 (I5.84p can be taken from Appendix O with 
Zk = Zi = z and the 7° and a 2 can be taken from Appendix |X1 In the nonrelativistic limit, 

z ^> 1, one gets 47 



64v^Fr 2V V 16 



46 It is the differential cross section for identical particles. The total cross section is J 4p2r 2 = inr 2 . 

47 This reproduces the result of Chapman and Enskog in the nonrelativistic theory for the shear viscosity. 
The vanishing value of the bulk viscosity is obtained in the limit m — > 00 |§j (Sections 8, 10). The result 
of the vanishing bulk viscosity of a monoatomic classical gas in the nonrelativistic theory is attributed to 



James Clerk Maxwell, see [77 1. 
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25 T 3/2 / 183 i * 
£ = = ^~ 3/2 1 + ... . (7.7) 



In the ultrarelativistic limit, z C 1, one gets 48 

3 T / I 



" = ^l 1 + ^ 2 + -l- < 78) 



? = 2^7^ { l + (I " 6 ln 2 + 67E ) z * + fc ' 111 2 + -) ' (7 9) 

where is the Euler's constant, 7# ~ 0.577. 

The perturbation of the distribution function (p (15. 5 7ft can be found too: 

V = na(T) + A ' T + AV ) V ^ + C°^%U^j , (7.10) 
where the C° is equal to 



= ^5_o(T) z 2 if 2 2 (^ 

" 64vr r 2 (15^ 2 + 2)K 2 (2z) + (3z 3 + 49z)fsr 3 (2z) 



(7.11) 



and the A is equal to 



, 2 _ 1 a(T)z 2 A- 2 2 (^)[(5- 37)^-37] 



64tt r 2 2fsT 2 (2z) + zK 3 (2z) ' K ' ' 

The A and the A 1 are used to satisfy the matching conditions (I5.65P and are equal to 

A o = A2 aV-(a 3 ) 2 ; Al = ^ oV-oV ^ ^ ^ alfl3 _ (q2)2 ^ (? 13) 

where the a s can be found in Appendix |A] In the nonrelativistic limit z>l one has 

* = z 2V - 2TW + 2ZT - + • ( 7 - 14 ) 

In the ultrarelativistic limit 2 C 1 one has 

V = ( -5*V + 8r - 12)V^ + 36^V, i [/J . (7.15) 

Note that although the shear viscosity diverges for T — > 00 the perturbative expansion over 
the gradients does not break down because the tp does not diverge (it tends to zero, con- 
versely) . 



48 The vanishing value of the bulk viscosity of a monoatomic classical gas in the ultrarelativistic limit is 
attributed to I. M. Khalatnikov, see 
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The phenomenological formula, coming from the momentum transfer considerations in 
the kinetic- molecular theory, for the shear viscosity is i] p h oc ln(\p\) (with the coefficient of 
proportionality of order 1), where (|p|) is the average relativistic momentum ( I4.2p . I is the 
mean free path. It gives the correct leading m and T parameter dependence of the (17. 5p 
with a quite precise coefficient 49 . The mean free path can be estimated as I ~ l/(a tot n) (see 
Appendix |D]). Choosing the coefficient of proportionality to match the nonrelativistic limit 
one gets [36( 

5 y/rtT K 5/2 (m/T) 
Vph 64v/? r 2 K 2 {m/T) ' 1 > 

If the bulk viscosity is expressed as £ p h oc ln{\p\) the coefficient of proportionality is not of 
order 1. In the nonrelativistic limit it is 25/(512v / 2z 2 ) and in the ultrarelativistic limit it 
is z 4 /(8647r). To reproduce these asymptotical dependencies the bulk viscosity should be 
proportional to the second power of the averaged product of the source term Q and the r 
that is to the (a 2 ) 2 . 

If a system has no charges, then terms proportional to the T\. in the (15.471) are absent. 
This results in very different values of a r k . In particular for the single-component gas in the 

case z ^> 1 one gets 

a 2 L l1 -n 4z 4 

l/LJ - +..., (7.17) 



a 2 | ?11 =i 15 



and in the case z <C 1 one gets 



2 I 

Q- |gii=o 

Q;2 1 9n=l 



13 + .... (7.18) 



In the case z ^> 1 this means for the bulk viscosity (with the constant cross section) that 
it becomes divergent over the temperature at least as 0{Tz 13 ^ 2 ) oc 0(T~ n / 2 ). However this 
dependence holds if the number- changing processes are not negligible. 

Although constant cross sections are the most simple and the most universal ones in 
approximate calculations, let's write down also formulas for some other simple energy de- 
pendencies of cross sections. Using the collision bracket from Appendix [Cj one gets for the 
cross section as /(2m) 2 = cru 2 /(2z) 2 (a is just a positive dimensional constant here) 

15T h 2 z 3 K 2 {z) 
m ~ 327ia 9z(3z 2 + 34)K 2 (2z) + (3z 4 + 157z 2 + 920)AT 3 (2;z) ' ^ 

c = 1 z^-3 1 )h-3 1 fK 2 {z) 

?1 32Tia{z 2 + 2$)Kz{2z) + QzK 2 {2 Z y 1 ' ' 

and for the cross section a(2m) 2 /s = a(2z) 2 /v 2 : 

„ _ 15T h 2 z 2 K 2 (z) 

V2 32vra (3^ 2 - 2)K 2 (2z) + z(3z 2 + l)K 3 (2z) ' 1 ' ' 



4!) 



This formula is justified only for rarified systems where the ideal gas equation of state is applicable. 
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32na (zK 3 (2z) - 2K 2 {2z)) 



(7.22) 



Using the low-energy current algebra isospin averaged tttt differential cross section [271 . 178 



O'CA 



3 647r 2 /* 



ml mi ( ml mi 
3 - 8^ + 7-f +11- 8^ + 16-£ 

s s 2 \ s s 2 



cos 2 e 



(7.23) 



{fn = 93 MeV is the pion decay constant) and treating all pions as identical particles at 
zero chemical potentials, one gets the formulas 



VCA 



360tt/4 



h 2 zK 2 (z) 



T 9z(93z 2 + 1730)K 2 (2z) + (69z 4 + 6167z 2 + 47104)AT 3 (2z) 



(7.24) 



P ^[(5-3 7 )fe-37] 2 ^) 

T (232 2 + 1024)iT 3 (2z) + 210zK 2 (2^)' 1 J 

which are exactly 2 times larger than the 1-st order calculations in the [27J 50 . Taking the 
scattering angle averaged cross section instead of the (I7.23p . one would get approximately 
the same result for the viscosities (errors are not more than 4%). In the low temperature 
limit, z>l, one gets the following expansions for the shear viscosities (I7.19p . (17.2 ip . (17. 24ft : 

*=^^( 1 -i I - 1+ -)- (726) 
" = ^i^( 1+ ^- 1+ -)- (727) 

120tt 3 / 2 fl , /2 / 2049 i \ . , 

and in the high temperature limit, 2<1, respectively: 

= J--^fl-^ 2 + -y (7-29) 



92tt a \ 23 



^ 2 = --^ 2 f 1 + ^(9 - A lE )z 2 - \z 2 \nz + ... ) . (7.30) 
7T a \ 2U 5 



45tt/*/ 17 j 



r? C A = — 1 z + ... . (7.31) 

92 T \ 368 



50 So that it looks like the double counting factor 1/2 is lost in the calculations of the viscosities in the 27 1 
(and presumably for the heat conductivity) , and is not lost in the current algebra total cross section. 
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In the low temperature limit, z> 1, one gets the following expansions for the bulk viscosities 

CL2QD, CL22D, UM- 

25 T q/9 / 247 i \ , 

6 = = -^ 3/2 1 z- 1 + ... , 7.32 

U 128^ a V 16 / 

25 T , /9 / 119 i \ . . 

' ~ 7 i--r^ +••• > ( 7 - 33 ) 



* 128v^Fa V 16 

^ G 23 T V 368 

and in the high temperature limit, z C 1, respectively: 

144U7T (T \ oU 

6 = t^-^ 2 ( 1 + ^(13 + 12 TB - 181n2> 2 + Az 2 lnz + ... ) , (7.36) 
367T a \ 6 

icA = ( X + 2l (89 + lUlE ~ UAln2 ^ 2 + 6z " lnz + -) • ( 7 - 37 ) 

Let's also write down the viscosities with the cross section er[l + c\v 2 j(2z) 2 + c 2 (2z) 2 /i; 2 ], 
unifying the previously considered ones: 

Vs = ^^Th 2 z 3 K 2 (z)/[z(2 + 306c 1 -2c 2 + 3(5 + 9c 1 + c 2 )z 2 )K 2 (2z) 

+ (920 Cl + (49 + 157 Cl + c 2 )z 2 + 3(1 + c x + c 2 )z 4 )K 3 (2z)}, (7.38) 

, _T ^[(5-3 7 )fr-3 7 ] 2 if 2 (^) 

^ 3 32vr ( r2(l + 3ci-c 2 )^ 2 (2z) + (20ci + (l + c 1 + c 2 )2 2 )ir3(22)' l ' ' 

where C\ and c 2 are some dimensionless coefficients such that c\ + c 2 + 1 > and ci > 
to make the cross section non-negative. The last formulas may have badly convergent 
expansions over the z and the z~ l for some values of the c\ and the c 2 so that they ain't 
expanded. 



B. The binary mixture 

The mixture of two species with masses mi, m 2 and the different classical elastic differen- 
tial constant cross sections af x , af 2 = c 2 i, 0" 22 is considered in this section. Using the ( 15.81)) 
with n 2 = and solving the matrix equation (I5.83j) one gets for the shear viscosity 

' 1 -[(xr^fC™, - 2x v x 2ll hlC% + (x 2/ j 2 ) 2 C™,], (7.40) 



10cr(T) A,, 
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where A,, = C%? V C°? 2 , - {C%f. The collision brackets for the C$ v fl5T85|) can be found in 
Appendix [C] and the 7° can be found in Appendix |A] 

In the important limiting case when one mass is large Z2 3> 1 (#2 and fa are finite so that 
Xi> <C 1) and another mass is finite one can perform asymptotic expansion of the special func- 
tions. Then one has X\i oc 0(1), 7° oc 0(1), x% oc 0(e~ Z2 z 2 ) 1 7° oc 0(z 2 )- The collisions 
of light and heavy particles dominate over the collisions of heavy and heavy particles in the 

C 2 ^ 2 z , and one has [ir^ir", 7r /J 7r I/ ] 2 i oc 0(z 2 ), C 2 ? 2 , oc 0(e * 2 2; 2 / ). In the Cy V the collisions of 

o 

light and light particles dominate, and one gets C®? v oc 0(1). And [jt^tt 17 , Wi^fu7]i2 cx 0(1), 
C®?2' oc 0(e~ Z2 z 2 ). In the shear viscosity the first nonvanishing contribution is the single- 
component shear viscosity ( 17. 51) . where one should take r 2 = cr^ and z = z\. The next 
correction is 



a 5 / 2 -*2 3Tg 2 e 

Ar 7 = z 2 e 2 — -=— — — — -. (7.41) 



64 v / 2^(3 + 32 1 + z 2 )^ 1 a 



The approximate formula 36l 



V = y^r/fca: fc , (7.42) 

where 77^ is given by the (17. 5p or the (I7.16P with mass and cross section c^, would give 
somewhat different heavy mass power dependence 0(e _Z2 zf ). 

Using the ( 15.921) with n\ = 1 and solving the matrix equation ( 15.911) one gets for the bulk 
viscosity 

p T (gya|) 2 T xi/x 2 /aia2 f74 ^ 
4 a(T) ^ a(T) " U " ' J 



Using the definition of the A\} 2 , (15.841) and the fact [r, ti]jw + [r, r]^ = (1C32j) one gets 



Ag/V = Xi'Xy [t, r]i2- Using the (I5.87P one gets [r, r] 12 > 0. Then using x^a\ + X2>oi\ = 0, 
coming from the (15.541) . the bulk viscosity can be rewritten as 

z = w , y 2 = j= v 1 > 0. 7.44 

The collision bracket [r, r] 12 can be found in Appendix [UJ and the a\ can be found in 
Appendix [A] 

In the limiting case z 2 3> 1 one has £1/ cx 0(1), £2' cx O^ -22 ^ 2 ), «i oc 0(e~ Z2 2:2 ), 
oc 0(1), A\\ cx cx 0(e~ Z2 Z2^ 2 ), [ r , r ]i2 oc 0(z2 1 ). Then for the bulk viscosity one gets 

= e -, 2/ /2 ff 2 Te-^ 2 ^ [2zf -5-2h* + 10/h] 2 

2 128v / 2^^icrf 2 (2i 2 + 3^ + 3)[2 2 -l-/i 2 + 5/ il ] 2 

VIII. CONCLUDING REMARKS 

The shear and the bulk viscosities of the hadron gas and the pion gas were calculated 
using the UrQMD cross sections. 
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The physics of the bulk viscosity is very interesting. In particular, in mixtures it can 
strongly depend on the mass spectrum. For instance, at the temperature T = 120 MeV 
(140 MeV) and zero chemical potentials the bulk viscosity of the hadron gas is larger in 
8.6-15.6 (14.6-40) times than the bulk viscosity of the pion gas. The used UrQMD cross 
sections have allowed to perform this comparison more accurately. 

Also the bulk viscosity can strongly depend on the quantum statistics corrections, the 
equation of state and the inelastic processes, which can be explained by nontrivial form of 
its source term(s). It's a future task to find the universal and optimal criterion for switching 
on/off the inelastic processes. Numerical calculations of the bulk viscosity along and around 
the chemical freeze-out line which don't involve the approximations of conserved or not 
conserved particle numbers (like in the 2ll . though the procedure of collisions of particles 



introduce some errors itself, which should be kept in mind) are desirable to get more accurate 
values of the bulk viscosity at these points. This is also needed for a better understanding of 
the chemical freeze-out itself and connected with it problems and to get a better description 
of the deviations from the chemical equilibrium. 

The transport coefficients are connected with fluctuations through the fluctuation- 
dissipation theorem. Because of strong dependence of the bulk viscosity on the equation 
of state one might expect to find its maximum at or near a phase transition point. Ac- 
cording to lattice calculations at zero chemical potentials [79| . the QCD phase transition is 
an analytical crossover. Calculations of the Polyakov loop |80j. the 't Hooft loop [8l| and 



some other calculations [82J suggest that hadronic degrees of freedom survive partially at 
some temperatures above the critical one. In the qualitative calculations for the hadron gas 
with the UrQMD total cross sections the maximum of the bulk viscosity has been found at 
T rs 190 MeV. This value fits in the transition temperature range T c = 185—195 MeV found 
from lattice calculations by the hotQCD collaboration, but calculations of the Wuppertal- 
Budapest collaboration with physical quark masses give the region T c = 150 — 170 MeV 
(83| . According to the Wuppertal-Budapest collaboration calculations of thermodynamic 
functions with physical quark masses [63|, the scaled trace of the energy-momentum tensor 
has its peak at T 190 MeV. In some other lattice calculations with somewhat different 
quark masses Him this peak is somewhat sharper and still is present at T « 190 MeV. 



According to the lattice calculations in the [63|, the squared speed of sound in the charge- 
neutral hadron gas (implying equal to zero and not developing chemical potentials) has it's 
minimum at T « 150 MeV, unlike the ideal gas calculations (T w 190 MeV) or the ones 



in the [19] (T 180 MeV). So that one could expect to find the maximum of the bulk 
viscosity somewhere in the temperature range T = 150 — 190 MeV, presumably closer to 
the lowest bound. This range covers the critical temperature values found from the lattice 
calculations. 

The shear viscosity is less dependent on the mass spectrum, the quantum statistics cor- 
rections, the equation of state and the inelastic processes. This may be explained by its 
more trivial source term(s). 
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Appendix A: The values of the a k , Y k and a' k 
Their definitions are 



a r k = (Q k , r r k ), Yk = frfc<<> n k ^k»), a% = (1, r fc s ) fc , (Al) 
where the round brackets are 

(F,G) k = 2 1 / F(p fe )C(p fe ) e -^. (A2) 



Then one can rewrite the 



l 



k zlK 2 {z k ) 



dr(r 2 - zl) l / 2 T s e- T . (A3) 



There is a recurrence relation for the a k \ 



ai = (s + l)a s - 1 + z^-s _ ( s _ 2 )^ar 3 . (A4) 
It can be derived from the (IA3[) written in the form 

< = -^j-t dr(r 2 - zlf^e- + z^. (A5) 
z k K 2 {z k ) J Zk 

Then after integration by parts the recurrence relation follow. Some values of the a k are 



1 



4 = -2(h k -4), (A6) 

z k 

4 = 1, (A7) 

a\ = h k - 1, (A8) 

a\ = 3h k + zl (A9) 

ai = (15 + z 2 k )h k + 2zl (A10) 

at = 6(15 + z 2 k )h k + z 2 k (l5 + zl), (All) 

at = (630 + 45*2 + 4)h + 5*2(21 + (A12) 
The a£ can be expressed through the a r k after the integration of the (I5.47[) (or the (I5.48j) if 



only the elastic collisions are considered) over the momentum, using the definition flAlj) . For 
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systems with only the elastic collisions some values of the a r k are written below, in agreement 
with g| (Chap. VI, Sec. 3): 

a° k = 0, (A13) 



, _ 2(c v -9)h k + 3hl-3z 2 k _ lk - 1 

a k — — 7> l Ai4 J 

C v 7fc ~ 1 



a\ = 2h k - 3^ - 3^^ = (5 - &y)h k - 3 7 fc^^, (A15) 

(+ 7fc 1 



where the assignments 7 and 7^ from [42] are used. They can be expressed through the c v 
and the c Vjk , defined in the (I5.44p . as 

7 = - + 1, 7fc = — + I- (A16) 

The 7^ can be rewritten as 

1 

Then it can be rewritten through the a r k : 



^ = li|7fS)/ ^-4f'^. (A17) 



Y t = ?(a' +4 - 2^a'+ 2 + (A18) 



Some values of the 7^ are 



7° = 10A fcj (A19) 
7^ = 10(6^ + ^), (A20) 
7^ = 10(7^ + ^(42 + ^)). (A21) 

Appendix B: The entropy density formula 

The Gibbs's potential is defined as 

$(P,T) = E(S,V) - ST + PV. (Bl) 
The differential of the energy is defined as 

dE = TdS - PdV + ^dN k = TdS - PdV + ^ fi a dN a , (B2) 
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where it is rewritten through the independent chemical potentials and the particle net 
charges N a . Then the differential of the $ reads: 

d$ = SdT + VdP + VadN a . (B3) 

a 

Because the $ is the function of the intrinsic variables P, T and the extrinsic N a the only 
possible form of it in the thermodynamic limit is 

$ = ^iV a a (P,T), (B4) 

a 

where 4> a are unknown functions. Then from the flB3j) one gets J^- = fi a , which means that 
4>a = Ha- Then substituting the (1B4|) into the (IBlj) one gets the relation 



N a fi a (P, T) = E{S, V)-ST + PV. (B5) 

a 

Being written for local infinitesimal volume it transforms into the expression 

KaHa = e-sT + P, (B6) 



from where the entropy density s can be found: 

e + P 



y^TWV (B7) 



T 

a 

Appendix C: The calculation of the collision brackets 

The momentum parametrization and the most transformations of the 12-dimensional 
integrals used below are taken from [42| (Chap. XI and XIII). Let's start from some assign- 
ments. The full momentum is 

p"=rt+tf l =rt+Ju = p" 1 - ( C1 ) 

The "relative" momentums before collision and after collision Q >fl are defined as 

<2 M = ^(Pku ~ Piw), = A£V fe „ - P'u,), (C2) 

with the assignment 

pH pv 

= <T - -pr, (C3) 

where P 2 = P^P^. The covariant cosine of the scattering angle can be expressed through 
the and the Q ,fl as 

cos e = Q ' Q '_ , (C4) 
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where ■ denotes convolution of 4-vectors. One also has Q 2 = Q' 2 and 



Q 2 = 4ml " (1 + «h) 2 P 2 = - (P 2 - M 2 kl ) 



P 2 V M; 



ki 



hi 



(C5) 



where 



= m k + mj, ^ w = 



m k + mi' 



Otkl = 



P 2 



sign(m fe - mi) J I 



M k , P 2 ' 



(C6) 



The function sign(x) is equal to 1, if x > and equal to —1, if x < 0. Note that not all P M 
and are independent: 



P»Q» = } POQ' =0. 



(C7) 



To come from the variables ip^iPii) ^° ^ ne variables (P M , Q M ) in the measure of integration 
first one has to come from the p^) to the (p^ + Pii,p k — Pi t ) (the determinant is equal to 
16) and then shift the relative momentum p k —Pxi on the a k iP^. Analogically for the (p' k , p'y) 
and the (P" J ',Q' fl ). The inverse relations for the P k ,Pii,p' k ,p'ii through the P^^Q^^Q" 1 are 



Pk 



±(i + au )p* + ±cy, 



(C8) 



p5i = ^(i-a«)j«'--g^, 



(C9) 



(CIO) 



'A 1 
Pli 



1(1 -aw)^-^. 



There is a need to calculate the following integrals 



(Cll) 



J, 



(a,b,d,e,f\q,r) 



Ikl 



ki 



TS(47ryz 2 z 2 K 2 (z k )K 2 ( Zl )a(T) 



-Q-Q' 

J^2 



kl- 



(C12) 



Let's start from the case of constant cross sections. After nontrivial transformations, de- 



scribed in more details in 42 , one arrives at 



J, 



(a,b,d,e,f\q,r) _ A d + e + WMkl 



(d,ej) q r [6/2] 



x 



X 



Zfc - gj 



zlzfK 2 (z k )K 2 ( Zl ) 

qi+r 1 +2k 3 

(-l) ri+fc2+fc3+h (2/i - 1)!! 



EE E E E^+ 2 < 

gi =0n=0 fc 2 =0 fc 3 =0 h=0 

b \ /„\ /„\ /rf+e 



,2(qi+ri+fc 2 +fc 3 ) 



(C13) 



/ 2(a + / - qt - n - & 2 - fc 3 ) + 3, 6 + 



d + e 



h+l,z k + zi 
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where 

d min(c(,e) 

*' /) = ^fr £ °™m*,g\ (ci4) 

where cr$ = 7^0^ is the classical elastic differential constant cross section. The is 
equal to the real, nonzero and non-diverging value (for any non-negative integer g) 

a < M = *i±± £ ixx , Pg{x) = {2g + i)___ZL___ i ( ci5) 

if the difference f — g is even and g < f. Above the -P g (x) is the Legendre polynomial. The 
K(d, e, g) is equal to the real, the nonzero and non-diverging quantity (for any non-negative 
integer g) 

d\ e \ 

K{< 3> = (d- 9 )!!(d + fl+ l)!i( e - fl )!l( e + 9+ l)IP ^ 

if g < min(d, e) and both the d — g and the e — g are even (which also implies that d + e is 
even). The [...] denotes the integer part. The integral / is 

/oo 
duu r+n K n (xu). (C17) 

Also there is the following frequently used combination of the J integrals 

j/(a,b,d,e,f\q,r) _ ^~~^^_^u ff \ ^2^2(/-u) j(a+fc,M+e,0,0|g+2u,r) _ j(a,b,d,e,f\q,r) (C18) 

The first term in the difference is obtained by the replacement of the Q' on the Q everywhere 
except for the W^i- Using this fact, the J' can be rewritten in the form 



, -I \ 1 1 (d,e,f) 1 r [6/2] 



?i=0ri=0 fc 2 =0 fc 3 =0 h=0 



x I V + 1 V f2(a + / - 9l - n - A; 2 - fc 3 ) + 3, 6 + - /i + 1, z k + z x I , 



where 



c l I min(d,e) 

a^ eJ) = ^(d + e + l)U(d + e,0,0)a^- £ K(d,e,g)a™ 

c l I min(d,e) 

= ^ I 1 - (d + e + 1) £ (d, e, ] . (C20) 
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There is a recurrence relation for the integral / (IC17I) [42] (Chap. XI, Sec. 1): 



I(r,n,x) = (r - l)(r + 2n - 1)1 (r - 2,n,x) + (r - l)x r+n - 1 K n (x) + x r+n K n+1 (x). (C21) 

For the calculations one needs only the integrals I(r,n,x) with the positive values of the 
n and the odd values of the r. If r > —2n + 1, the I integrals can be expressed through 
the Bessel functions K n (x), using the (1C21I) , when r = 1 or r = —2n + 1. Then using the 
recurrence relation for the K n (x) [86| 

2n 

K n+1 (x) = K n _ x {x) + —K n (x), (C22) 

the final result can be expressed through a couple of Bessel functions. If r < — 2n — 1, then 
the recurrence relation ( 1C21I) becomes singular if one tries to express the I(r,n,x) through 
the I(—2n + l,n, x). Using the flU2ip . the / integrals with r < — 2n — 1 can be expressed 
through the integrals G(n,x) 

/oo 
duu~ n - x K n {xu). (C23) 

There is a recurrence relation for the G(n,x): 

G(n,x) = -— (G(n- l,x) -x~ n K n {x)). (C24) 
2n 

It can be easily proved by the integration by parts of the flC23j) and using the following 
relation for the K n (x) j86| 

o n 

—K n (x) = —K n (x) - K n ^(x). (C25) 

OX X 

It is found that collision brackets have the simplest form if they are expressed through 
G(n, x) with n = 3 or n = 2 and the Bessel functions K 3 (x) and K 2 (x) or K 2 (x) and Ki(x). 
It was chosen to take G(x) = G(3, x) and K 3 (x), K 2 (x). The G(x) can be expressed through 
the Meijer function [87J 



G(x) = ^Gl'° 3 ( {x/2f 



1 

-3,0,0 



(C26) 

The needed scalar collision brackets can be expressed through the J' as 



' V s V^/' 1 N \«/ r \f S I T/(0,r+s-M-w,M,u,0|r-it,s-u) /'r'OV^ 
u=l v=l • * 



r r si 1 \ % \ ^ / r \ / S \ T /{0,r+s-u-v,u,v,0\r+s-u-v,0) fnoQ\ 

u=l w=l V / \ / 
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and the needed tensorial collision brackets can be expressed as 

[rW,^« — (* )( T n )(- — 2 .0|2 + n 1)2+ n 2 ) 

ni=0n 2 =0 V 2/ V 1/ 
_|_ 2 j/(l,«i+«2,r-ni,s-n2,l|l+«i) 1 +"2) _|_ j7(0,ni+n2,r-ni,s-7i2,2|ni,7i2)\ 

1 A s + A + 1 V i \s+l-n 2 / T/(l,ni+n 2 ,r+l-ni,s+l-n2,0|l+ni,l+n 2 ) 

ni=0n 2 =0 V * / \ 1 / 
T/(0,7ii+n 2 ,r+l-Tii,s+l-Ti2,l|rii,n2)\ , 2r r+2 *+2i , 1 2r r+2 si 

+ \4V r ,tI +2 \u - ^zlzf[r r ,rl] kl , (C29) 

[rW, r^] kl ~EE( n lt) (/ (2,n 1+ n 2 ,-„ llS -„ 2 ,0|4 + „ 1+n2 ,0) 

ni=0n 2 =0 ^ ' ^ ' 
— 2j'( 1 ' ni+n2,r ~ ni,S ~ n2,1 ^ 2+ni+n2 ' ^ + j'(°: n i+ n 2,'"-rii,s-ri 2 ,2|rii+n 2 ,0)\ 



r+l s+1 

/(l,ni+n2,r+l— ni,s+l— n 2 ,0|2+ni+n 2 ,0) 

2^+3 - m n2 ; ^ ni j v w 

ni=0n 2 =0 ' x 

T/(0,m+n 2 ,r+l-7ii,s+l-Ti2,l|m+ri2,0)\ . 2r _ +2 s +2t , ^ Jr r+2 _sl 
~ J ki ' 3*- * k 3 ' T J fc ' 

+ ^V^-^VV (C30) 

Below some lowest orders collision brackets are presented with the following notations: 

~ K 3 (z k + zi) ~ _ K 2 {z k + zi) ~ G(z k + zi) 



K 2 (z k )K 2 ( Zi y K 2 {z k )K 2 { Zi y K 2 (z k )K 2 ( Zi y 

Z k i = z k + 2/, = ^fe — ^- (C31) 

For the scalar collision brackets one has: 

-It, nU = [r, r] kl = °f'* (Pg^i + P^K 2 + P™K 3 ), (C32) 

where 

P^ = -2Z„(4 + 44^ - 2Z* ), (C33) 

PS' 1 ) = 4(3Z 2 + 8) + 324^ + 8Z£ h (C34) 

= -HA (C35) 



and 

r 21 _ r_2 _1 _. ff fc* ^ (p^K 4- P {1 ' 2) ^ 4 P (1 ' 2) /< 1 

r^ijfc^-r ' ^ fc ~ 4z 2 z 2 z 2 ^ g11 312 s13 ' 

V / k I kl 



(C36) 
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P£? = 4i(Zti + 102% + 16) + 12z 5 kl Z kl (Z 2 kl - 24) - ^(Z* - 72Z 2 kl - 384) 
- 192z 3 kl Z kl (Z 2 kl + 40) - 484^(13^, - 224) + 96*^(7^, - 80) 
+ 48^(13^ + 48), (C50) 

= -z^Z^zUzl - 6) + 12z kl Z kl - 6Z 2 kl }. (C51) 
And for the tensor collision brackets one has: 

[W 7 , = °L 79 / y272 (P^?Ki + P^K 2 + P$S>K,), (C52) 
a{i ) nz k z l /j kl 



where 



and 



P^ } = -2Z kl [z e kl (5Z 2 kl - 8) + 24afc(Z* - 16) - lUz 2 kl Z 2 kl (Z 2 kl + 8) 

+ 48Z*(Z* +72)], (C53) 

P(° 12 0) = 4(52j - 40Z* - 64) - 244(5ZJ + SZ 2 kl + 128) 

+ 576z 2 kl Z 2 kl (Z 2 kl - 16) - 192^,(52* + 16), (C54) 

4°i3 0) = -5«[4(^ - 24) - 24^], (C55) 



^]« = 4l 79 Z 272 (P&^i + 4°22 0) ^ + Pff^s), (C56) 



where 



Pff = 2Z H [4(8 - 5Z 2 kl ) + 72z 4 kl (3Z 2 kl - 8) - A80z 3 kl Z kl (Z 2 kl - 4) 

— 33§z kl Z kl (Z kl + 8) + 240^^ + 8) + 192Z£,(Z& + 67)], (C57) 

PS 2 0) = 4(5Z^ - A0Z 2 kl - 64) + 2404^ 3 / - 244(5Z* + A8Z 2 kl - 192) 

+ 1920z 3 kl Z kl (Z 2 kl " 8) - 1924^(17^ - 112) + 1920 z kl Z 3 kl (Z 2 kl - 8) 

+ 768^4(5^ + 6), (C58) 

4 23 0) = -54^[4(^/ - 24) + A8z kl Z kl - 24Z fc 2 J. (C59) 

If z k = zi, then the G(x) function is eliminated everywhere and the collision brackets simplify 
considerably. 

The constant cross sections are not the only possible ones resulting in the analytical 
expressions for the collision brackets and, hence, for the transport coefficients. The analytical 
expressions through the Bessel and the Meijer functions (possibly, relatively simple ones) 
can be obtained also using some non-constant cross sections. For example, this is possible in 
the case when cross sections are proportional to integer powers of the P 2 and/or the cosO. 
As one can see from the ( IC12j) . an integer power of the P 2 would result just in shift of the 
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index a which can be easily taken into account in the ( IC13j) or the (1C19j) . A power of the 



cos would result in shift of the first index of the a^' 9 ^ (1C15|) . and one would get 



min(d,e) 

"'■'-'l''> -H±L- \ (f + l)- 1 - (d + e + i) K(d,e,g)<rV + ™\, (C60) 

9=0 



a(T) 



instead of the olf' (1C20j) . So that the J or the J' integrals just get a different factor. 



Half-integer powers 51 of the P 2 can also be taken into account in the same way as the 
integer powers of the P 2 however one would need to introduce one more special function 
through which to express the I integrals with even r using the recurrence relation (|C21j) . 
The convenient choice of this function is found to be 

/CO 
duu- n+2 K n (xu). (C61) 

Then one can express these functions through the function (the n = 3 is found to be the 
convenient choice) 



G 2 (x)^G 2 (3,x) = ^-- 4 Glj[(x/2) 2 



1 

_3 3 q 
2 ' 2 ' 



(C62) 



using the recurrence relation 



G 2 (n,x) = ^—^[K n (x)x~ n+Z - G 2 (n - l,x)}. (C63) 

Its derivation is the same as the derivation of the ( 1C24I) . Using the above mentioned pre- 
scriptions one can calculate collision brackets for a quite large class of cross sections. Few 
examples of simple energy dependencies for cross sections are considered below. 

Using the cross sections auiP 2 / (T 2 Z| ; ) = auiv 2 jZ\ h growing with the energy, one finds for 
the scalar collision brackets {pu is just a positive dimensional constant here, and = jki&ki) 

-[r,n] kl = [r,r]„ = ^-^(P^K 1 + P^K 2 ), (C64) 

°V ) z k z i 



where 



and 



^ii M|1) = \z kl (4i - HA + Zfa + 48 ^)> ( C65 ) 



Ps2 ' ^ — ~ z ti ~ z ki^ki + 2Z^, (C66) 



[r-rflu : [""'■ ~i ]//, -^^^(P^ ll) K 1 + P^K 2 , (COT) 



<T) z\z 2 Zl x 



51 Actually all powers of the P 2 can also be taken into account in the same way if all mutual differences 
between all powers of all the P 2 terms are integers. 
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where 



Pi!i 2|1) = Z kl (4i + HiZm + - 104^ - 5*h3h + 8 ^(^ + 30)), (C68) 



Pll? m = ^[-4(^ + 8) + 4^(^-16) + 24^(^-16) 

- 24^(^ + 8) - ^,(^2, + 8) + ZUZI + 80)], (C69) 



and 



r_ . 1 ... 4?r /p(i,2|i)g- . p(i,2|i)y , 

where 

)(1,2|1) _ -7 /_,,5 , r> 4 7 /|J 7 2 _ 1f . 2 7 3 , r„. v4 , o 7 3/ 7 2 



and 



where 



and 



where 



(C70) 



} = Zki(-4i + HiZki - 44^ - 104^ + 5^ + 8Z* (Z£ + 30)), (C71) 



^22 2|1) = \\4i{Zl + 8) + 4^(^-16) -24^(^-16) 

- 2z 2 kl ZUZ 2 kl + 8) + z H Z* (2* + 8) + Z*(Z* + 80)], (C72) 



[J J, _ °ff ^ / p (2,2|l)~ p(2,2|l)y p (2,2|l)yN , p7 oX 

[T ,r!jfct - ~~^7f\ Q z 2 z 2 Z 4 i^sll K l + ^sl2 K 2 + ^13 A 3j, ^7dj 
V / k I kl 



P? n m = 2[-zt l Z kl (5Z 2 kl + 4) + 18z 4 kl Zl l (Z 2 kl -2) 

- 18z 2 kl Z 5 kl (Z 2 kl + 70) + 6Z* (Z* + 222Z 2 kl + 5120)], (C74) 

^ s ( i 2 2 2|1) = 4(-^ + 44^ + 32) + 364^(3^ + 8) 

- 364^(11^ - 8) + 12Z° (19Z* + 856), (C75) 

P$? W = 4A (C76) 



r 2 2l _ °"fei ^ /p(2,2|l)y , p(2,2|l)^ p (2,2|l) ^ x ^--v 

V / k I kl 



P s ( 2 2 i 2|1) = 2Z fcI [4(5Z& + 4) + 124Z„(^ - 10) - Qz^ZUZl - 46) 
- 2±z 3 kl Z 3 kl (Z 2 kl + 20) - 64^(^ + 182) + 12z u Z^ + 98) 
+ 6Z* (3* + 226^ + 5440)], (C78) 
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P&? 1] = 4i(Zti - 44^ - 32) + 2Azt l Z kl {Zl l + 40) + 12^(5^ - 184) 
- 48«(7^ - 80) - UzlZ^Zl + 184) + 24z„Z* (13Z£ + 40) 
+ 12^(19^ + 904), (C79) 

Pg£ ll) = ~4A (C80) 
And for the tensor collision brackets one has: 

f^.^Im = -^*> { p»*m> + i«f>& + if* >*,>, (csi) 

where 

^ti'i'^ = ^ z ti z ki(Zli - 1) + 9^jZ|j(Z^ - 2) 

- 9*^(2* + 22) + 3Z*(Z* + 1263* + 2240)], (C82) 



4i 2 0|1) = 4(-5^ - 142* + 16) - 184^(3^ - 8) 

- 18z 2 kl Z 4 kl (5Z 2 kl - 8) + 6ZUl3Z 2 kl + 376), (C83) 



4°i3 0|1) = HAn (C84) 



and 



p=3Z=37 - ^ 2?r /p(0,0 1 y , p(0,0|l)^ , p (PfiWfr\ (CRZ,\ 



where 



P4 0|1) = 2Z kl [-2zt l {Zl - 1) + 154Z„(Z£ - 4) + 64^(^2 + 23) 
- 304^(Z^ + 8) - HiZWzl + 614) + 15z H Z* (2* + 68) 
+ 6Z* (2Z* + 377^ + 8480)], (C86) 

4 22 0|1) = 4(5^ + 142* - 16) - 60z 5 kl Z kl (Z 2 kl - 8) + e^Z* (29Z* - 184) 
- 2404^(^ - 8) - 6z 2 kl Z 4 kl (S5Z 2 kl + 184) + 60z kl Z 5 kl (5Z 2 kl + 8) 
+ 6^(67^ + 2824)), (C87) 

p£f> = -54^- (C88) 
Using the descending cross sections OkiZlJv 2 one finds for the scalar collision brackets 



_[ T , Tl]a = [ T , T]u = Jl * (pWI-^ + fA'l-)^ + pi»'->/f 3 ), (C89) 

<j{i ) oz k z l 
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where 

= 2Z kl (-4i + 84 + 8^), (C90) 

= 4i(Z 2 kl ~ 8) - 84(3^ + 8) - 64Z 2 kl , (C91) 

^g' lhl) = 4^("4(^ - 24) + 24Z* ), (C92) 

and 

r 21 _ r 2 i _ ^ 71 ( p(l,2|-l) £^ , p (l,2|-l) ^ p(l,2|-l) ^ \ /p Q o\ 

F,TlJfcI — LT > Tljlfc - ~^(rj 64~zV ^ 511 Al + ^12 ^2 + ^13 ^3j, (^93) 

where 

^ii 1 ! 2 '"^ = ^(-24^ + 324 + 6Az 3 kl Z kl + 5124 + 256z kl Z kl + 512^), (C94) 

P^' l) = z 5 kl Z kl (Z 2 kl - 16) - 164(^ + 8) - S2z 3 kl Z kl (Z 2 kl + 8) 

- 2564(^4 + 8) - 128z kl Z kl (Z 2 kl + 8) + 128Z fc 2 ^ - 16), (C95) 

P^- 1] = z^ZU-zUZl - 32) + 16z„Z„ + 32Z*), (C96) 

and 

[r,r 2 ] H = [tV] h = ^^(P^K, + P™~»K 2 + P™~»K 3 ), (C97) 

where 

Pgf- l) = 2Z kl {z\ l Z kl + 164 - 324Z H + 2564 - l28z kl Z kl + 256Z 2 kl ), (C98) 

PSa hl) = 4iZki(Z 2 kl - 16) - 164(^ + 8) + S2zl l Z kl (Z 2 kl + 8) 

- 256z 2 kl (Z 2 kl + 8) + 128z kl Z kl (Z 2 kl + 8) + 128Z 2 kl (Z 2 kl - 16), (C99) 

P£s 2hl) = 4iZU4i(Z 2 kl - 32) + 16z kl Z kl - 32Z 2 kl ), (C100) 



and 

fj T 2i„_ ^ ^ / P (2,2|-i)y . p (2,2|-i)y p (2,2|-i)yx mnn 

L - - 3840^2^2 l-^sll A l + -^12 -^2 + ^13 A 3j, (U1U1J 

where 

p ( 2 ,2|-i) = 2Z«[4(7^- 76) - 604(5^, -32) 

- 115204^ - 6) + S840Z 2 kl (Z 2 kl + 18)], (C102) 
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PS? hl) = 4i(-7Z 4 kl + 188Z 2 kl + 320) + 604(5Z fc 4 z + 80^ - 256) 

+ 230404(Z 2 - 24) + 7680Z 2 (7Z 2 - 72), (C103) 

P& 2hl) = 4^(4(7^4 - 300) - 30022,), (C104) 

and 

\J1 21 _ °kl 71 ( p(2,2|-l) g- p (2,2|-l)^ (2,2|-1)^ v m 0^ 

F ,T \kl - ^7^ 3840^2 ^21 ^1 + ^22 K 2 + ^ s 23 K 3>, (L>IU&) 

where 

P s ( 2 2 i 2hl) = 2Z„[-4(7Z£ - 76) - 3604Z« + 604(5Z 2 + 64) - 11520z 3 kl Z kl 

- 38404(^ - 22) + 76S0z kl Z kl (Z 2 kl - 10) + 3840Z 2 (Z 2 kl + 22)], (C106) 

^ ( 22 2M) = 4(7^ - 188^ - 320) + 1204Z„(3^ + 16) 

- 604(5^ + 48Z 2 + 512) - 38404^(^4 - 24) 

- 76804(7^ + 88) + 15360z kl Z kl (Z 2 kl + 40) + 7680^(7^ - 88), (C107) 

P££ hl) = 4^[4(300Z H -7^)- 1204(3^ -32) 

+ 60^Z H (5^- 32) + 3840Z2j. ( C108 ) 

And for the tensor collision brackets one has: 

= -^^-^(^fl-'lif, + + PS^'K,), (C109) 

where 

P^f 1 ^ = 2Z kl [zU7Z 2 kl - 76) - 604(5^ - 32) - 115204(Z 2 - 6) 

+ 3840^(^ + 18)], (C110) 

4i2° hl) = 4i(-7Zt t + 188^ + 320) + 604(5^ + 80^ - 256) 

+ 230404(Z 2 - 24) + 7680Z 2 (7Z 2 kl - 72), (Clll) 



4°i3° hl) = 4^(4(7^ - 300) - 3002%), (C112) 



and 



1^=37 - °A 71 /p(o,o|-i)^ , p(Q,o|-i)y i p(Q,Q|-i)g-\ rrn^ 

[TT^TT , TV^Im - ^j^^o^V T21 Al + / T22 A 2 + -f+23 A 3j, (Ulld) 

where 

= 2Z fe ,[-4(7Z 2 - 76) -18004Z fcJ + 604(5Z 2 - 32) 

- 38404(7Z 2 Z - 22) + 19200zhZ h (Z* - 4) + 76802* (22* + 11)], (C114) 
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4 22 hl) = 4i(7Zt t - 188Z 2 kl - 320) + 600z 5 kl Z kl (3Z 2 kl - 16) 

- 604(5^ - - 256) - 384004^ + 76804(17Z* - 88) 

- 76800 z kl Z kl {Z 2 kl - 8) + 7680Z^(13Z^ - 88), (C115) 

P^~ l) = 4^[4(300Z w -7^)-6004(3Z fc V64) 

+ 300^^ - 64) + 38400Z*]. (C116) 

Also averaged cross sections with the considered energy dependencies can be written down 
easily. Below there are only averaged energy insertions: 

, * /7 2 S _ r7 - 2 4?' ' '°' 0|0 ' 0) = Mz4 + z l + + ^K 2 {z kl )) 

{ / kl) kl JWAOICO) Zkl{Z 2 M _ z l l)K ,{Z kl ) + ^lK 2 {Z kl ) ' ^ n 



(Z 2 kl /v 2 ) = Z, 



J, 



(-1,0,0,0,0|0,0) 
kl 



kl ^(0,0,0,0,010,0) 
J kl 

U&ZhGjZu) + Z kl (Z 2 kl - z 2 kl )K 3 (Z kl ) - 4(4 + Z 2 kl )K 2 (Z kl ) 

Z kl (Z 2 kl -z 2 kl )K 3 (Z kl )+Az 2 kl K 2 (Z kl ) ■ 1 ) 

Using these averaged insertions one could reproduce correct asymptotic temperature de- 
pendencies of the viscosities however with a somewhat different coefficient in the high- 
temperature limit. 

For complicated forms of cross sections one can calculate numerically the 1- or 2- 
dimensional integrals (corresponding to the unknown energy and angular dependencies). 
As long as implemented UrQMD cross sections don't depend on the scattering angle (being 
treated as averaged ones over it) one has to calculate only 1-dimensional integrals numeri- 
cally. For this purpose the J integrals can be represented in the form 



[6/2] 

7 2(a+f)+b+(d+e)/2-h+5 
J kl 



T (a,b,d,e,f\q,r) _ 7l(d + e + ^ h f b \ 2 

J « - ztz?K 2 (z k )K 2 ( Zl ) ^-^-Ur* 

/oo 
du*? k f\z klU )(u 2 + Zkl /z kl y{u 2 - z kl /z kl yu 2 ^-«-^ b - h -^ 2 

An 2 - ir^' 2 +f+\u 2 - 4/^) (d+e)/2+/+1 ^ + e)/ 2+ 5-, + i(^), (C119) 
where the a[ d k f' f \z k iu) is generalized to 

d , . mm(d,e) 

^i J \Zkiu) = Yl ° {f ' 9) me,9), at\(Z kl u) = lkl a kl (Z kl u). (C120) 

And the J' integrals can be represented in the form 

[6/2] 



T ,(a,b,d,e,f\q,r) _ 7r(c? + e - 1)!! >^ fc / 6 \ 2 (o+/)+6+(*fe)/2-h+5 

/oo 
^^(z^)^ 2 + z kl /z kl y{u 2 - z kl /z kl yu 2 ^-^ +b - h -^i 2 

An 2 ~ l^W+V - 4/^) (d+e)/2+/+1 ^ +e) /2 + 5-, + i(^), (C121) 
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where the a kl ' e (Z k iu) is generalized to 

d( 7 \ ( min(d,e) \ 

°tf J \Zuu) = I 1 - (d + e + 1) g jr(d, e, j . (C122) 

To calculate collision brackets faster one can bring all integrated expressions under one 
integral and simplify the integrand. 



Appendix D: The collision rates and the mean free paths 

The quantity o k % >o d 3 p' 'ud 3 p' \», which enters in the elastic collision integral (15.221) . 

p k ,p ii ,p k /P if 

represents the probability of scattering per unit time times unit volume for two parti- 
cles which had momentums p^ and pu> before scattering and momentums in the ranges 
(p'khP'k' +dp' k ,) and {p'n^p'ui + dp' iv ) after the scattering. The quantity gk'j^jsfk' rep- 
resents the number of particles per unit volume, which have momentums in the range 
(Pk' , Pk' + dpk>) ■ The number of collisions of particles of the k'-th species with particles 
of the I'-th species per unit time per unit volume is then 52 

Zel ^ nn ill' f d 3 Pk' d 3 PlV d 3 P ' k , d 3 p' u , (o) (0) , , 

aw = gygi'jK-xs / — 5 — io — _,o .0 h' hv w uv- l m J 

\^) J P k ' PlV P k' P XV 

To get the corresponding number of collisions of particles of the A;'-th species with particles 
of the Z'-th species per unit time per particle of the k'-th species, Rfi V , one has to divide 
the (ID1I) on the ^kn'^k' (recall that oc gy by definition), which is the number of particles 
of the /c'-th species per unit volume divided on the number of particles of the k'-Xkx species 
taking part in the given type of reaction (2 for binary elastic collisions, if particles are 
identical, and 1 otherwise). This rate can be directly obtained averaging the collision rate 
with fixed momentum pk of the fc-th particle species 

rr,ei _ f d 3 Pw ,3 , , 3 / ,(o) W k v , , 

' = m ' lk " J Wjs' 1 p « d MM' m 

over the momentum with the probability distribution j^f^ (and spin states which is triv- 
ial): 



w _ .. .. Ik'v f d 3 Pk ' d 3 p iv d 3 p' k , d 3 p' iv (Q) (0) R% 



K k'v = 9k'9i'j—^ — / — o o lo ~h' J XV w wi' • [ U ' i 

{2-K) b n k > J Pl> Pl v p' k , p'\ v im'ny 

So that to get the mean rate of the elastic collisions per particle of the fc'-th species with 
all particles in the system one can just integrate the sum of the gain terms in the collision 
integral (I5.22p over , d )3 Pk — and average it over spin: 

i$ = l> fe V (D4) 



It represents some sum over all possible collisions. In the case of the same species one factor jk'i' just 
cancels the double counting in momentum states after scattering and another factor ^k'v also reflects the 
fact that scattering takes place for ("*') ss \r?y pairs of undistinguishable particles in a given unit volume. 
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One can express the R kn , through the j^ ' < ' ' l >°) integrals from Appendix [C] as 

Del /rp\ ,-(0,0,0,0,010,0) m r-\ 

For simplicity let's consider constant cross sections in what follows. Then the ( 1D5[) becomes 
Rt'v = 9k'9i'lk'l' — H — [(**' - z v ) 2 K 2 {z k t + z v ) + z k ,z v {z k > + z v )K^(z k > + z v )\, (D6) 

71° 

where a k , v is the classical elastic differential constant cross section of scattering of a particle 
of the k'-th species on particles of the Z'-th species. For the case of large temperature or 
when both masses are small, z k < "C 1 and Zj( < 1, one has expansion 

~ Arr cl T 6 / 1 \ 

Rtv = gwai'-rwv-^- 1 1 - i (4 + 4) + ■■■) • (D7) 

For the case of small temperature or when both masses are large, z k > ^> 1 and Zy ^> 1, one 
has expansion 

Sei „ y/2a% v T*z k ,z V y/zy + z^e~ z k'- z i' ( %z\ + 19z k ,zy + 8zp 



Rtn> = m 5/2 1 + I 7 1 7 + - • (D8) 

7T 5 / 2 V Sz k /Zy{z k ' + Zv) ) 

For the case when only one mass is large, zy 3> 1, one has somewhat different expansion 

~ e , y^grT^ + l)*?/V«*'-«" / 44 + 15^+15 \ 

*W - <7*«* ^ + h> ■■■) ■ (D9) 

The a(T)4^°'°' 0|0 ' 0) in the « can be replaced in the limit of high temperatures 
with Anaf/^dvk'l) and in the limit of low temperatures with 4:ira kll , ( | vy \ ) \A + m k > ~]~my = 
47rcTfc/j'(|£VeZ,fc'Z'|)j where {\v k /\) is the mean modulus of particle's velocity of the k'-th species, 

Jd 3 Pk^\p k >) 2e -«*'(l + * fc ,) mim 

\l^'l) = — — — mt, — — = -Trrrr-\ - \ — — ts~f~l \ ' l mu ) 



Jd 3 p k >fl°\p k >) 4> K 2( z k') \ nz k , K 2 (z k 

and (|tVei,A;'Z'i) is the mean modulus of the relative velocity, which coincides with the (|#fc'|) 
for high temperatures. Then the resultant collision rate Rf, would reproduce simple non- 
relativistic collision rates know in the kinetic-molecular theory. To get the (approximate) 
mean free time one has just to invert the R e k 1 ,: 

4 = 4- (DH) 
K k' 

The (approximate) mean free path if, can be obtained after multiplication of it on the 

l'i = { W- (D12) 

n k> 

For the single-component gas one gets 

jei = Ml ire-^(z 1 + l) 

Lv Rf, v g^af^zlK^zxY [ 6) 
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The nonrelativistic limit of the (1D13|) with the gi = 1 coincides with the same limit of the 
formula 



If = <!*!> = I (D14) 

ATTaf ini (\v re i\) foaftmyfi' 

which is the mean free path formula coming from the nonrelativistic kinetic-molecular theory 
obtained by Maxwell. The ultrarelativistic limit of the (lD13j) with the g\ = 1 coincides with 
the same limit of the formula 

Analogically one can introduce inelastic rates R™ el of any inelastic processes to occur for 
the k'-th particles species. Then the mean free time t l $ el in what any inelastic process occurs 
for the particles of the k'-th. species can be introduced as 



The mean free path for the particles of the fc'-th species is obtained through the rate 
Rf, + R l y el and can be written as 



Dei i TDinel ' 

(D17) 
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